ESC-TR-95-096 


Quarterly  Technical  Report 


Solid  State  Research 


1 9960621  094  i995:4 


Lincoln  Laboratory 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
Lexingtois,  Massachusetts 


Prepared  for  the  Department  of  the  Air  Force  under  Contract  F19628-95-C-0002. 


Approved  for  public  release;  distribution  is  unlimited. 


DSC  01 


This  report  is  based  on  studies  performed  at  Lincoln  Laboratory,  a  center  for 
research  operated  by  Massachusetts  Institute  of  Technology.  Thework  was  sponsore<l 
l)y  the  Department  of  the  Air  Force  under  Contract  F19628-95-C-0002. 

This  report  may  be  reproduced  t<j  satisfy  needs  of  U.S.  Government  agencies. 


The  ESC  Public  Affairs  Office  has  reviewed  this  re[)ort,  and 
it  is  releasable  to  the  National  Technical  Information  Service, 
where  it  will  be  available  to  the  general  j>ublic,  including 
foreign  nationals. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication, 
FOR  THE  COMMANDER 


Contracted  Support  Management 


Non-Lincoln  Recipients 

PLEASE  DO  NOT  RETURN 

Permission  is  given  to  destroy  this  document 
when  it  is  no  longer  needed. 


MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
LINCOLN  LABORATORY 

SOLID  STATE  RESEARCH 

QUARTERLY  TECHNICAL  REPORT 


1  AUGUST-31  OCTOBER  1995 

ISSUED  17  JUIVE  1996 

Approved  for  public  release;  distribution  is  unlimited. 


LEXINGTON 


MASSACHUSETTS 


ABSTRACT 


This  report  covers  in  detail  the  research  work  of  the  Solid  State  Division  at  Lincoln 
Laboratory  for  the  period  1  August  through  31  October  1995.  The  topics  covered  are 
Electrooptical  Devices,  Quantum  Electronics,  Materials  Research,  Submicrometer  Technol¬ 
ogy,  High  Speed  Electronics,  Microelectronics,  Analog  Device  Technology,  and  Digital 
Integrated  Circuits.  Funding  is  provided  primarily  by  the  Air  Force,  with  additional  support 
provided  by  the  Army,  ARPA,  Navy,  BMDO,  NASA,  and  NIST. 
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INTRODUCTION 


1.  ELECTROOPTICAL  DEVICES 

A  simple  scheme  has  been  developed  for  calculation  of  the  Strehl  ratio  and  mode  coupling  effi¬ 
ciency.  Accurate  large-numerical-aperture  microlenses  and  precision  optical  alignment  are  needed  for 
efficient  coupling  of  diode  lasers  to  single-mode  fibers  or  waveguides,  and  this  scheme  is  useful  for  quick 
diagnosis  of  the  optical  system  and  helpful  in  guiding  the  microlens  design  and  fabrication. 

2.  QUANTUM  ELECTRONICS 

First,  second,  third,  fourth,  and  fifth  harmonics  of  the  output  of  a  fiber-pumped  passively 
g-switched  Nd:YAG  microchip  laser  have  been  obtained  at  pulse  energies  of  8.0,  3.5,  0.3,  0.7,  and 
0.01  jXS,  respectively,  in  an  optical  head  occupying  a  volume  <  3  cm^.  This  compact,  economic,  all-solid- 
state  source  provides  coherent,  subnanosecond,  multikilowatt  pulses  in  the  infrared,  visible,  and/or  ultra¬ 
violet  at  repetition  rates  in  excess  of  10  kHz. 

3.  MATERIALS  RESEARCH 

High-quality  Pbj_^Uj^Te/PbTe  multiple  quantum  wells  (MQWs)  have  been  grown  by  molecular 
beam  epitaxy,  and  their  measured  Seebeck  coefficients  at  300  K  were  found  to  be  about  twice  that  of  the 
best  bulk  PbTe.  This  is  the  first  detailed  experimental  investigation  of  MQW  structures  designed  to 
improve  the  thermoelectic  figure  of  merit  ZT,  and  the  value  of  >1.2  achieved  for  these  PbTe 

quantum  wells  has  broken  the  decades-long  ZT  =  I  barrier  for  a  room-temperature  thermoelectric 
material. 

4.  SUBMICROMETER  TECHNOLOGY 

Plasma-deposited  carbon-based  polymer  films  have  been  developed  for  use  as  a  dry-deposited,  dry- 
developed  positive-tone  photoresist  for  193-nm  lithography.  These  films  crosslink  upon  exposure  to 
193-nm  radiation,  enabling  selective  silicon  uptake  via  reaction  of  silylamine  gas  with  hydroxyl  in  the 
film. 


Field  emission  of  electrons  has  been  demonstrated  from  nitrogen-doped  diamond  at  fields  <0.1 
V  jum“^  Although  the  vacuum  field  is  small,  these  electrons  have  energies  of  several  kilovolts  and  are 
emitted  from  the  surface  in  well-defined  beamlets. 

5.  HIGH  SPEED  ELECTRONICS 

GaAs  complementary  metal-insulator-semiconductor  field-effect  transistors  (MISFETs)  with  a  low- 
temperature-grown  GaAs  gate  insulator  have  been  fabricated  using  the  same  epitaxial  layer  structure  with 
self-aligned  Si  and  Be  implants.  The  maximum  drain  current  is  40  and  120  mA/mm  for  a  normally-off 
n-  and  p-channel  MISFET,  respectively,  and  increases  to  500  mA/mm  for  a  normally-on  «-channel 
device. 
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6.  MICROELECTRONICS 


The  signal  variance  of  photoelectrons  collected  in  the  well  of  an  electronically  shuttered,  back- 
illuminated  charge-coupled  device  has  been  measured  to  be  substantially  less  than  the  expected  shot  noise 
of  the  incident  light.  Under  certain  conditions  there  is  a  factor  of  1.75  improvement  in  the  squared  signal- 
to-noise  ratio. 

7.  ANALOG  DEVICE  TECHNOLOGY 

Third-order  intermodulation  power  generated  by  a  Josephson  junction  in  a  YBa2Cuj0.7_^  micro- 
wave  resonator  has  been  measured  and  compared  to  a  calculation  based  on  the  resistively  shunted 
junction  model.  The  results,  which  are  relevant  to  the  coupled-grain  model  used  to  describe  microwave 
properties  of  high-T^  superconductors  and  to  the  understanding  of  Josephson  junctions  in  these  materials, 
agree  with  junction  parameters  determined  by  power-dependent  loss  measurements  on  the  same 
resonator. 

8.  DIGITAL  INTEGRATED  CIRCUITS 

A  decoded  version  of  the  ultradense,  once-programmable,  read-only  memory  has  been  designed 
and  is  in  fabrication.  The  program  goal  is  to  build  a  memory  with  bit  dimensions  approaching  the 
lithographic  pitch,  and  components  such  as  single  bits  at  0.25-/tm  diameter  and  undecoded  arrays  have 
been  demonstrated. 
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XXIV 


1.  ELECTROOPTICAL  DEVICES 


1.1  SIMPLE  SCHEME  FOR  ANALYSIS  OF  DIODE  LASER  MICROOPTICS 


The  interest  in  microoptics  for  use  with  diode  lasers  has  been  increasing.  Microoptics  can  be  used 
for  collimating  the  laser  output  or  for  efficient  coupling  of  optical  powers  into  single-mode  fibers  or 
waveguides  [1].  Because  of  the  relatively  large  divergence  of  the  output  from  a  diode  laser,  large- 
numerical-aperture  microlenses  are  needed  and  high  degrees  of  accuracy  are  also  required  for  the  optical 
alignment  [1].  In  this  work,  we  use  a  simple  analytical  technique  that  is  capable  of  quick  diagnosis  of 
the  optical  system  and  is  highly  desirable  for  the  microlens  design  and  fabrication. 

Figure  1-1  shows  a  collimating  microlens  in  front  of  a  point  source.  By  equating  the  optical  path 
lengths  SP  and  SOQ,  one  can  readily  derive  the  lens  profile: 
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where  n  is  the  refractive  index  of  the  lens  material.  This  is  the  well-known  hyperbolic  profile  for  exact 
collimation  of  a  point  source  without  any  aberration.  At  large  numerical  apertures,  however,  small 
deviations  from  this  ideal  lens  profile  or  inaccuracies  in  the  alignment  to  the  source  can  result  in  severe 
aberrations.  In  this  work,  we  first  consider  a  special  type  of  deviation  that  is  important  for  mass-trans¬ 
ported  microlenses  [2]--[4],  but  the  analysis  technique  is  rather  generally  applicable  to  other  types  of 
errors. 
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Figure  1  -L  Microlens  for  collimation  of  a  point  source.  By  considering  the  optical  path  lengths^  phase  distortion 
in  the  plane  PP '  can  readily  be  calculated  for  any  deviation  from  the  ideal  lens  profile  and  the  precise  optical 
alignment.  The  phase  deviations  can  in  turn  be  used  to  predict  the  lens  performance. 
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One  frequent  error  in  the  mass-transported  microlenses  is  the  size  of  the  lens  sag,  which  arises  from 
inaccurate  etched  mesa  height  in  the  fabrication  of  the  lens  preform  [2], [3].  When  the  preform  is  overetched 
by  a  small  fraction  e,  the  lens  profile  becomes 


z'  =  (l  +  e) 
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where  is  the  designed  focal  length.  For  small  numerical  apertures,  this  simply  results  in  a  change  of 
effective  focal  length  from  /q  to 

r  =  /„/(l+e)  .  (I'S) 


without  significant  degradation  in  lens  performance.  However,  if  the  source  is  placed  at  the  correct  focal 
length/,  the  output  is  only  approximately  collimated.  For  larger  numerical  apertures,  there  is  an  increas¬ 
ing  phase-front  distortion,  which  can  be  calculated  by  evaluating  the  optical  path  length  difference 

5  =  s-f--nz'  ,  (1-4) 

with 


^  =  ^(f'+zf+y^  . 

The  calculated  phase  variation  across  the  lens  aperture  can  in  turn  be  used  to  predict  the  lens  performance. 
The  solid  trace  in  Figure  1-2  shows  the  calculated  phase  variation  for/g  =  160  /tm  and  e  =  0.1.  Note  that 
the  phase  deviation  has  a  predominantly  y‘*  dependence,  as  can  be  seen  in  the  logarithmic  plot  in 
Figure  1-3. 

To  assess  the  impact  on  lens  performance,  consider  a  cylindrical  lens  for  simplicity.  One  measure 
of  the  effect  of  aberrations  of  a  lens  is  the  Strehl  ratio  SR,  which  is  the  ratio  of  the  peak  far-field 
intensities  with  and  without  phase-front  distortions.  For  uniform  illumination  of  a  cylindrical  lens,  the 
Strehl  ratio  is  given  by 


where  D  is  the  lens  aperture  and  A  is  the  wavelength  in  air.  The  Strehl  ratio  is  an  indicator  of  the 
efficiency  with  which  light  is  coupled  into  subsequent  stages  of  an  optical  system,  e.g.,  coupling  into  a 
fiber. 

Equation  (1.6)  can  be  evaluated  as  a  function  of  the  numerical  aperture  with  e  as  a  parameter,  as 
illustrated  in  the  solid  traces  in  Figure  1-4.  Note  that  a  Strehl  ratio  close  to  1  is  maintained  for  e  =  0.02, 
but  severe  degradation  results  at  large  numerical  apertures  when  e=  0.1.  It  should  also  be  noted,  however, 
that  the  Strehl  ratio  for  e  =  0.1  can  be  considerably  improved  by  fine  adjustment  of  the  source  position 
for  the  introduction  of  an  opposite  wavefront  curvature,  as  illustrated  by  the  dashed  traces  in 
Figures  1-2  and  1-4.  However,  freedom  for  such  adjustment  may  not  always  be  available  in  some 
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y(Mm) 


Figure  1-2.  Phase  distortions  of  mass-transport  microlens  that  has  been  overetched  by  10%  in  the  fabrication  of 
the  lens  preform.  The  parameter  fg  is  the  designed  focal  length  of  the  microlens  without  the  overetch. 


in  y 


Figure  1-3.  Logarithmic  plot  of  phase  distortion  as  a  function  of  radial  position  resulting  from  10%  overetch  of 
the  lens  preform  (cf.  the  solid  trace  in  Figure  1-2).  This  plot  clearly  shows  a  predominantly  y"^  dependence  of  the 
phase  deviation. 
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Figure  1-4.  Calculated  Strehl  ratio  in  the  case  of  uniform  illumination  for  cylindrical  mass-transport  microlenses 
with  two  different  overetches  in  the  preform  fabrication. 


systems,  such  as  anamorphic  microlenses  for  astigmatism  correction.  Further,  note  that  the  Strehl  ratio 
for  a  circular  lens  is  more  degraded  (relative  to  a  cylindrical  lens)  because  the  area  weighting  in  the 
integral  analogous  to  Equation  (1.6)  emphasizes  the  more  distorted  outer  regions  of  the  circular  lens. 
Although  this  can  also  be  analyzed,  the  following  discussion  will  still  be  limited  to  cylindrical  lenses  for 
consistency. 

For  nonuniform  lens  illumination,  the  coupling  efficiency  between  more  general  mode  profiles  can 
readily  be  shown  to  be 
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where  is  the  amplitude  profile  of  the  incident  beam  after  lens  collimation,  ytg  is  the  amplitude  profile 
of  the  mode  to  be  coupled  into,  and  N  is  the  appropriate  normalization  factor.  (Note  that,  in  actual 
coupling  of  a  laser  mode  to  a  single  mode  fiber  or  waveguide,  an  additional  focusing  lens  or  its  equivalent 
is  needed.  The  mode  profiles  in  Equation  (1.7)  are  those  at  the  plane  of  the  collimating  lens.)  By 
comparison  with  Equation  (1.7),  the  Strehl  ratio  given  by  Equation  (1.6)  can  be  interpreted  as  the 
coupling  efficiency  between  modes  of  uniform  intensities.  Conversely,  Equation  (1.7)  can  be  considered 
as  a  weighted  Strehl  ratio  for  a  nonuniformly  illuminated  lens. 
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For  a  case  of  practical  interest,  we  now  consider  the  coupling  between  two  Gaussian  modes  of  the 
same  width,  i.e., 


T]  = 


J-D/2 


(1.8) 


where  is  the  amplitude  profile  of  the  incident  beam.  Two  rather  typical  diode  laser  output  beam 
divergences  of  38  and  28°  full  width  at  half-intensity  maximum,  with  corresponding  a’s  of  90  and  66  fim, 
respectively,  are  used  for  the  calculation  shown  in  Figure  1-5.  The  coupling  efficiencies  start  at  low 
values  for  small  numerical  apertures  because  the  incident  light  is  not  fully  captured.  As  the  numerical 
aperture  increases,  the  efficiencies  increase  and  become  higher  than  those  of  the  uniformly  illuminated 
lenses  because  of  less  contribution  from  the  more  distorted  outer  regions. 
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Figure  1-5.  Calculated  coupling  ejficiency  in  the  case  of  a  beam  with  Gaussian  profile  in  one  dimension  incident 
on  cylindrical  mass-transport  microlenses  with  two  different  preform  overetches  and  two  different  incident  beam 
divergences. 


The  present  technique  is  rather  generally  applicable  to  other  types  of  aberrations  and  is,  in  fact,  not 
limited  to  microlenses.  As  an  example,  we  have  carried  out  calculations  for  aberrations  resulting  from 
small  misalignment  in  the  y-direction  (i.e.,  with  the  source  slightly  off-axis)  of  a  uniformly  illuminated 
cylindrical  lens,  as  shown  in  Figure  1-6.  Note  that  the  problem  is  more  severe  for  lenses  made  of  lower- 
index  material.  It  is  also  worth  noting  that  tolerances  greatly  improve  by  using  the  alternate  lens  configu¬ 
ration  [5]  as  shown  in  Figure  1-7. 

Z.  L.  Liau 
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Figure  1-6  .  Calculated  Strehl  ratios  for  slight  misalignment  of  the  source  off  the  optical  axis  of  uniformly 
illuminated  cylindrical  lens.  The  calculation  has  been  carried  out  for  two  different  misalignments  and  for  lenses  of 
different  refractive  indices. 
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Figure  1-7.  Alternative  lens  configuration  that  offers  superior  tolerances  for  lens  and  alignment  accuracies. 
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2.  QUANTUM  ELECTRONICS 


2.1  PASSIVELY  e-SWITCHED  MICROCHIP  LASERS  FOR  ULTRAVIOLET 
GENERATION 

Experimental  studies  of  passively  <2-switched  picosecond  microchip  lasers  [1]  and  the  fourth- 
harmonic  conversion  of  the  output  of  these  devices  into  the  ultraviolet  (UV)  have  been  described  previ¬ 
ously  [2].  Since  these  reports,  we  have  generated  the  third  and  fifth  harmonic  of  the  microchip  laser 
output,  performed  further  characterization  of  the  device,  and  worked  out  the  theory  of  harmonic  fre¬ 
quency  conversion.  These  advances  are  discussed  here,  following  a  brief  review  of  the  previously  re¬ 
ported  work. 

The  principle  behind  the  operation  of  a  passively  Q-switched  laser  is  that  the  presence  of  a  satu¬ 
rable  absorber  within  the  laser  cavity  prevents  the  onset  of  lasing  until  the  average  inversion  density 
within  the  cavity  reaches  a  critical  value.  The  onset  of  lasing,  at  that  point,  produces  a  high  intracavity 
optical  field  which  quickly  bleaches  the  saturable  absorber,  increasing  the  Q  of  the  cavity  and  resulting 
in  a  g-switched  output  pulse.  The  short  cavity  length  of  a  microchip  laser  results  in  short  optical  pulses, 
and  an  infrared  passively  Q-switched  Nd:YAG/Cr:YAG  microchip  laser,  pumped  by  a  fiber-coupled 
1-W  CW  diode  laser,  has  produced  pulses  as  short  as  218  ps  (full  width  at  half-maximum  [FWHM])  with 
peak  powers  up  to  25  kW  and  ideal  beam  quality  (diffraction-limited  TEMqq  transverse  mode,  transform- 
limited  single  longitudinal  mode,  linearly  polarized).  It  is  the  high  peak  power  and  ideal  beam  quality  of 
the  passively  (2-switched  infi'ared  microchip  laser  that  allows  for  efficient  frequency  conversion  to  the 
visible  and  UV  regions  of  the  spectrum  using  nonlinear  optical  crystals. 

The  0-switched  microchip  laser  used  for  the  harmonic  generation  discussed  below  had  a  total 
cavity  length  of  1.25  mm,  comprising  1  mm  of  Nd^''':YAG  and  0.25  mm  of  Cr‘^''':YAG,  and  an  output 
coupler  with  85%  reflectivity.  It  generated  8-/J,  440-ps  pulses  at  a  repetition  rate  of  10  kHz.  The  output 
at  1.064  ftm  was  single  frequency  and  diffraction  limited,  with  a  6.8-mrad  single-sided  divergence  and 
a  l/e^  waist  diameter  of  100  urn.  Although  we  have  constructed  passively  ^-switched  microchip  lasers 
with  higher  peak  powers,  this  device  had  a  combination  of  peak  power  and  waist  size  (and  hence  confocal 
parameter)  that  were  well  matched  to  the  requirements  for  efficient  UV  generation  in  the  desired  compact 
format  without  any  intermediate  optics. 

The  1.064-^m  radiation  from  this  microchip  laser  was  frequency  doubled  with  44%  energy  effi¬ 
ciency  in  a  5-mm-long  crystal  of  KTP  placed  adjacent  to  the  laser  (with  no  intervening  optics)  and 
oriented  for  type  n  phase  matching.  The  532-nm  output  had  a  circularly  symmetric,  nearly  Gaussian 
(super-Gaussian)  profile  with  a  5.9-mrad  divergence  and  an  85-fim  waist  diameter.  The  3.5-fiJ  green 
pulses  had  a  duration  of  380  ps. 

To  generate  the  UV  harmonics  of  the  microchip  laser  output,  we  used  antireflection-coated  pieces 
of  P-B3B2O4  (BBO)  oriented  for  type  I  phase  matching.  In  all  cases,  the  crystals  were  5  mm  long,  had 
flat  faces,  and  were  placed  in  the  output  path  of  the  laser  with  no  intervening  optics.  Before  discussing 
the  experimental  UV  generation,  we  will  briefly  review  the  theory  of  harmonic  generation  in  BBO,  in 
the  limit  of  low  conversion  efficiency. 
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BBO  has  significant  double  refraction.  For  third-,  fourth-,  or  fifth-harmonic  generation  of  the 
output  of  the  Nd:YAG  microchip  laser,  the  harmonically  generated  beam  walks  out  of  the  aperture 
defined  by  the  generating  beams  in  a  length  that  is  short  compared  to  the  length  /  of  BBO  used.  The 
aperture  length  for  coherent  interaction  of  the  beams  is  [3] 

7 

^  tan(p3)  ’  (2.1) 


where 


1  1 
-  +  - 


(Oel 


-1/2 


(2.2) 


in  which  w  is  the  1/e^  radius  of  a  nearly  Gaussian  beam,  the  subscripts  1 ,  2,  and  3  refer  to  the  generating 
and  generated  beams,  respectively,  and  the  subscript  e  indicates  that  the  waist  dimensions  are  measured 
in  the  direction  perpendicular  to  the  ordinary  optic  axis  of  the  crystal  (o  indicates  the  direction  parallel 
to  the  ordinary  axis);  and  where  the  walk-off  angle  Pg  is  given  by 


tan(p3) 


”o3  1 

9  ^ 
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.«e3  . 
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in  which  and  are  the  ordinary  and  extraordinary  refractive  indices  at  the  frequency  of  the  gen¬ 
erated  beam  and  6  is  the  phase-matching  angle  in  the  crystal  defined  by 


M  =  2d 


Aj  ^2  A] 


=  0  , 


(2.4) 


with  A  being  the  ffee-space  wavelength  and 

n-1/2 
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The  peak  intensity  of  the  generated  beam  is  given  by 


(2.5) 


(2.6) 
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where 


^eff  =  <^3 1  ® + P3  )  “  ^22  ^ + P3  )  sin  3^ 


(2.7) 


is  the  effective  nonlinear  coefficient  (0  =  90°  for  a  properly  oriented  crystal),  c  is  the  speed  of  light,  e 
is  the  dielectric  constant  of  free  space,  and  JiS6)  <  1  is  a  function  of  the  amount  of  power  in  one  of  the 
generating  beams  that  falls  within  the  acceptance  angle  of  the  crystal.  When  the  acceptance  angle  of  the 
crystal 


ft  — 
t'acc  ~ 


2k 


dAk 


1-1 


(2.8) 


is  smaller  than  the  divergence  angle  6^  (in  the  BBO)  of  a  generating  beam,  we  will  use 
^  ^acc 

(2.9) 

For  all  calculations,  we  will  assume  IdgjI  «  \d22\,  1^22!  =  2.2  pnW  [4],  and  the  refractive  indiees  given 
in  Table  2-1  [5]. 


TABLE  2-1 

Refractive  Indices  of  /3-BaB204  at  Nd:YAG  Harmonics 


X 

(nm) 

"0 

1064 

1 .54254 

1.65510 

532 

1 .55552 

1 .67493 

355 

1 .57757 

1 .70556 

266 

1.61461 

1 .75707 

213 

1.67467 

1 .84707 

The  peak  power  and  pulse  energy  £3  in  the  generated  beam  are  calculated  from  the  peak 
intensity.  In  the  dimension  parallel  to  the  ordinary  axis  of  the  BBO,  we  will  assume  that  the  generated 
beam  has  a  nearly  Gaussian  profile  with  a  l/e^  radius  of  In  the  orthogonal  dimension,  the  direction 
of  walk-off,  the  profile  is  roughly  top  hat,  with  a  diameter 


^3  - 


/ 

tan(p3)  ■ 


(2.10) 
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The  peak  power  is  approximately 

V2/3 

’  (2.11) 

and,  after  accounting  for  the  temporal  shape  of  the  pulse,  the  pulse  energy  is  approximated  by  [6] 

£3=1.16P3t3  ,  (2-12) 

where  a  reasonable  estimate  of  the  pulse  duration  (FWHM)  is 
-1/2 

(2.13) 

The  green  output  of  the  KTP  was  doubled  to  produce  266-nm  UV  radiation  with  20%  energy 
efficiency  in  a  5-mm-long  crystal  of  BBO  placed  adjacent  to  the  KTP.  The  far-field  profile  of  the 
resulting  0.7-/J  UV  pulses  was  nearly  Gaussian  in  the  dimension  parallel  to  the  ordinary  axis  of  the  BBO, 
with  a  l/e^  single-sided  divergence  of  4.1  mrad.  In  the  opposite  dimension  the  far-field  profile  approxi¬ 
mated  a  sin^x/x^  pattern,  with  an  angular  separation  of  1.1  mrad  between  the  two  central  nodes.  The  near¬ 
field  and  temporal  profiles  of  the  UV  output  were  not  measured. 

From  the  theory  presented  above,  for  fourth-harmonic  generation  of  the  microchip  output  6  =  47.3°, 
P3  =  -4.84°,  =  1.62  pnW,  0)^3  =  30  pm,  =  423  pm,  =  269  ps,  =  627  pm,  =  1.687, 

and  J{5e^)  =/(502)  =  0.91.  With  /j  =  4  =  287  MW/cm^  (/^  =  this  gives  calculated  values  of 

4  =  17.9  MW/cm^,  P3  =  2844  W,  and  £3=  0.89  pJ;  in  excellent  agreement  with  the  0.7  pJ  obtained.  From 
the  theoretical  calculations,  the  far-field  divergence  of  a  diffraction-limited  beam  in  the  direction  parallel 
to  the  ordinary  axis  would  be  2.8  mrad,  indicating  that  the  actual  beam  was  ~1.4  times  diffraction  limited, 
similar  to  the  532-nm  generating  beam.  In  the  opposite  dimension,  the  far-field  profile  of  a  423-pm-diam 
top  hat  should  have  central  nulls  spaced  by  1.2  mrad,  in  good  agreement  with  the  observed  value. 

To  generate  the  third  harmonic  of  the  1.064-pm  microchip  output,  the  fundamental  frequency  was 
combined  with  the  second  harmonic  in  a  properly  oriented  BBO  crystal  placed  adjacent  to  the  KTP. 
When  the  KTP  was  oriented  for  maximum  532-nm  generation,  the  output  energy  at  355  nm  was  0.24  pJ/ 
pulse.  The  output  was  increased  to  0.3  pJ/pulse  by  rotating  the  KTP  in  the  plane  normal  to  the  propa¬ 
gation  of  the  1.064-pm  radiation.  The  far-field  profile  of  the  355-nm  output  (under  conditions  of  maxi¬ 
mum  output  energy)  was  qualitatively  similar  to  that  obtained  at  266  nm,  with  a  1/e^  divergence  of 
3  mrad  in  the  nearly  Gaussian  dimension,  and  a  central  node  spacing  of  1.0  mrad  in  the  orthogonal 
dimension. 

When  the  KTP  was  oriented  for  maximum  532-nm  generation,  the  peak  intensity  of  the  infrared 
remaining  after  the  doubling  crystal  was  112  MW/cm^  [7(1.064  pm)  —  7(532  nm)].  However,  as  a  result 
of  the  type  II  phase  matching  in  the  KTP,  only  half  of  this  was  properly  polarized  to  contribute  to  third- 
harmonic  generation  in  type  I  phase-matched  BBO,  with  the  result  that  7j  =  56  MW/cm^  and  =  287 
MW/cm^.  Under  these  conditions,  the  calculated  pulse  energy  at  355  nm  is  0.26  pJ  (to  be  compared  with 


h  - 


.2  ■'■.2 
h  h 
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the  measured  value  of  0.24  fjJ).  By  properly  rotating  the  KTP,  the  intensities  of  the  first-  and  second- 
harmonic  radiation  are  brought  into  balance,  while  the  fraction  of  the  first  harmonic  that  is  properly 
polarized  for  third-harmonic  generation  is  increased.  As  a  result,  the  amount  of  third  harmonic  generated 
was  increased  by  25%. 

The  266-nm  radiation  was  also  combined  with  the  fundamental  to  generate  213-nm  light  in  a 
properly  oriented  crystal  of  BBO  placed  adjacent  to  the  frequency-quadrupling  crystal.  When  aligned  for 
maximum  output  energy,  the  ordinary  axis  of  the  fifth-harmonic  crystal  was  perpendicular  to  the  ordinary 
axis  of  the  fourth-harmonic  crystal,  and  fifth-harmonic  pulse  energies  of  ~  0.01  /J  were  obtained. 
(Because  of  the  low  average  power  at  this  wavelength,  the  uncertainty  in  the  energy  measurement  is 
~  50%.)  The  far-field  profile  of  the  213-nm  output  was  qualitatively  similar  to  that  of  the  266-nm  light, 
except  rotated  by  90°  due  to  the  relative  orientations  of  the  fourth-  and  fifth-harmonic  crystals.  The 
l/e^  divergence  in  the  nearly  Gaussian  dimension  was  ~  1.5  mrad,  and  the  central  node  spacing  in  the 
other  dimension  was  0.8  mrad. 

As  in  the  case  of  third-harmonic  generation,  only  half  of  fundamental  remaining  after  second- 
harmonic  generation  in  the  KTP  is  properly  polarized  to  contribute  to  fifth-harmonic  generation.  In  this 
case,  the  useful  1.064-;Um  radiation  is  polarized  perpendicular  to  the  ordinary  axis  of  the  BBO  crystal 
used  for  fourth-harmonic  generation,  and  experiences  double  refraction  in  that  crystal,  with  a  “walk-ofP’ 
angle  of  ^.03°.  As  a  result,  it  overlaps  the  fourth-harmonic  radiation  at  a  position  of  uniform,  maximum 
intensity.  If  we  ignore  the  divergence  of  the  1.064-jUm  radiation,  the  theory  presented  above  predicts  a 
pulse  energy  of  0.01  /xJ.  At  the  position  of  the  fifth-harmonic  crystal,  the  fundamental  beam  should  have 
diverged  by  40%  in  radius,  decreasing  the  amount  of  generated  fifth  harmonic  proportionally. 

We  also  built  a  single  device  that  simultaneously  generated  the  first  five  harmonics  of  Nd:YAG. 
This  was  done  by  positioning  a  third-harmonic  crystal  immediately  after  the  KTP,  after  the  fourth- 
harmonic  crystal,  or  after  the  fifth-harmonic  crystal.  Different  positions  of  the  third-harmonic  crystal 
resulted  in  different  output  powers  at  each  of  the  UV  wavelengths. 

The  most  thoroughly  tested  of  the  devices  discussed  above  is  a  frequency-quadrupled  passively 
g-switched  microchip  laser.  The  laser  and  the  nonlinear  crystals  are  packaged  in  a  stainless-steel  cylinder 
1  cm  in  diameter  by  3  cm  in  length  [2].  As  a  result  of  the  short  aperture  length  (IJ  in  the  BBO,  the 
temperature  window  for  UV  generation  is  large  and  at  least  1  mW  of  time-averaged  UV  output  is 
obtained  over  a  temperature  range  from  0  to  70°C.  Currently,  this  device  has  been  operated  for  >  10^*^ 
shots  without  degradation. 

The  use  of  5-mm-long  crystals  of  BBO  in  all  of  the  experiments  described  above  is  somewhat 
arbitrary.  When  the  aperture  length  in  the  nonlinear  crystal  is  short  compared  to  the  crystal  length,  as  is 
the  case  here,  the  conversion  efficiency  is  linearly  dependent  on  the  crystal  length,  to  the  extent  that  pump 
depletion  and  beam  divergence  can  be  neglected.  The  amount  of  UV  generated  at  any  of  the  Nd:YAG 
harmonics  could  be  increased  simply  by  using  a  longer  piece  of  BBO.  Since  the  acceptance  angle  and 
temperature  window  for  nonlinear  conversion  are  determined  by  the  aperture  length,  they  are  not  affected. 

J.  J.  Zayhowski  C.  Cook 
C.  Dill  ni  J.  L.  Daneu 

J.  Ochoa 
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3.  MATERIALS  RESEARCH 


3.1  fflGH  THERMOELECTRIC  FIGURES  OF  MERIT  IN  PbTe  QUANTUM  WELLS 

Recent  renewed  efforts  to  improve  thermoelectric  cooling  materials  have  stimulated  considerable 
interest  in  the  concept  of  creating  multiple-quantum-well  (MQW)  superlattices  to  increase  thermoelectric 
properties  [1].  Calculations  by  Hicks  and  Dresselhaus  [2]  showed  a  very  large  increase  in  the  two- 
dimensional  (2D)  thermoelectric  figure  of  merit  ^2^  of  a  single  Bi2Te3  quantum  well  as  the  thickness 
of  the  well  decreased,  because  of  the  large  density  of  states  per  unit  volume  that  occurs  for  small  well 
widths.  Here,  we  report  a  first-of-its-kind  investigation  of  the  thermoelectric  properties  of  Pbj_^u^Te/ 
PbTe  MQW  structures  grown  by  molecular  beam  epitaxy  (MBE)  on  BaF2  substrates,  where  enhanced 
thermoelectric  properties  are  also  anticipated  via  the  same  effect.  The  PbEu-chalcogenide  material  system 
is  well  suited  for  investigating  thermoelectric  properties  of  quantum  wells,  because  it  is  both  a  good 
thermoelectric  material  [3]  and,  unlike  Bi2Te3,  has  a  good  epitaxial  growth  technology  base  [4]-[6]. 

The  figure  of  merit  for  a  thermoelectric  material  is  given  by 
ZT  =  ^(friK  =  PpTiK  ,  (3.1) 

where  S,  o,  K,  T,  and  Pp.  (=  S^d)  are  the  Seebeck  coefficient,  the  electrical  conductivity,  the  thermal 
conductivity,  the  temperature,  and  the  thermoelectric  power  factor,  respectively.  High  Pp  and  low  K  are 
desirable  to  obtain  the  highest  possible  ZT  for  efficient  cooling  and  energy  conversion  applications.  For 
homogeneous  material.  Equation  (3.1)  can  also  be  written  as 

ZT  =  S^neiJ.TI(Kj+  K^)  ,  (3.2) 

where  n,  e,  fi,  X),  and  are  the  carrier  density,  electronic  charge,  carrier  mobility,  and  the  lattice  and 
electronic  parts  of  the  thermal  conductivity,  respectively.  The  electronic  part  of  the  thermal  conductivity 
in  Equation  (3.2)  is  determined  from  the  Wiedemann-Franz  law  as 

K  =  ’  (3-3) 

where  is  a  parameter  dependent  on  the  electron  scattering  mechanism  and  the  statistics  [7].  Here,  we 
used  =  4.4  X  IQ-^  Q  W/K.  For  PbTe  we  used  ^  =  20  mW/cm  K  [8]. 

In  the  case  of  an  MQW  stmcture,  if  we  assume  that  current  and  heat  flow  only  through  the 
quantum  wells,  we  get  a  2D  thermoelectric  figure  of  merit 

Z^^T  =  S\eilTI(icr  +  C)  •  (3-4) 

where  n^,  Kp',  and  are  the  carrier  density  and  lattice  and  electronic  parts  of  the  thermal  conductivity 
of  the  quantum  wells,  respectively.  The  2D  power  factor  is  likewise  given  by 
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(3.5) 


where  tr^  is  the  electrical  conductivity  of  the  quantum  wells. 

The  MBE  growth  was  carried  out  on  1.8  x  1.8-cm,  as-received,  cleaved  BaFj  (111)  substrates  with 
thicknesses  of  ~  1  mm,  which  were  thermally  cleaned  in  the  MBE  system  at  490°C  for  32  min  prior  to 
growth  of  a  Pbo  95gEuo  042X0  buffer  layer  at  350°C.  Reflection  high-energy  electron  diffraction  (RHEED) 
patterns  indicated  a  change  in  growth  morphology  from  3D  islands  to  2D  layer-by-layer  growth  and  a 
smooth  surface  for  subsequent  MQW  structure  growth.  Buffer-layer  growth  was  then  continued  an  ad¬ 
ditional  8  min  for  a  total  thickness  in  the  range  1100-1900  A.  Since  the  accuracy  of  the  measured 
thermoelectric  properties  improves  with  film  thickness,  relatively  thick  MQW  structures  in  the  range 
4.5-6  jjm  were  grown  at  about  300°C  with  the  number  of  periods  typically  equal  to  120.  Growth  rates 
were  in  the  range  0.75-1.0  /im/h.  Effusion  cells  containing  Eu,  zone-refined  PbTe  [9],  and  Te2  were  used 
in  growing  the  undoped  PbQ  95gEuQ  042'!'^  buffer  layers.  Only  one  effusion  cell  containing  zone-refined 
PbTe  was  used  for  the  quantum-well  growth,  but  four  effusion  cells  were  used  for  the  Bi-doped 
Pbo927EuQ023Te  barriers,  i.e.,  PbTe,  Eu,  Te2,  and  Bi.  Good-quality  RHEED  streaks  were  observed  through¬ 
out  the  growth  of  the  MQWs  indicating  high-quality  2D  growth  for  the  entire  structure. 

A  schematic  cross  section  of  the  Pbg  927Euo.o73^®^^^^®  MQW  structure  along  with  the 
Pbo95gEuoo42Te  buffer  layer  and  the  BaF2  substrate  is  shown  in  Figure  3-1.  Thin,  undoped  buffer  layers 
were  used  so  that  the  electrical  conductance  of  the  buffer  layer  would  be  negligible  relative  to  the 
conductance  of  the  approximately  120  quantum  wells  in  each  film.  The  electrical  conductance  and 
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Figure  3-L  Schematic  cross  section  of  the  PbEu-chalcogenide  multiple-quantum-well  structure  investigated. 
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Seebeck  coefficients  of  the  MQW  stractures  were  measured  with  the  temperature  gradient  and  electric 
field  applied  parallel  to  the  epilayer  interfaces.  Theoretical  calculations  by  Hicks  and  Dresselhaus  [10] 
indicate  that  the  electrical  conduction  through  the  barrier  layers  is  negligible  compared  to  that  through 
the  quantum-well  layers  for  the  MQW  structures  used  in  the  experiments  reported  here.  Also,  the  con¬ 
ductance  of  the  thin,  low-doped,  low-carrier-mobility  Pb^  ^Eu^Te  buffer  layer  was  negligible.  Numerous 
Pbj_^u^Te/PbTe  MQW  samples  were  grown  with  PbTe  well  widths  in  the  range  10-50  A  and  Pbj  ^u^Te 
barrier  widths  in  the  range  300-600  A.  Samples  with  x  =  0.042  exhibited  a  large  decrease  in  electron  mobility 
with  increasing  barrier  width,  indicative  of  3D  conduction.  For  the  case  of  26-A  wells  and  1300-A  barriers, 
the  300-K  carrier  mobility  for  a  carrier  density  of  ~  10^^  cm“^  (averaged  over  the  well  and  barrier)  was 
only  240  cm^A^  s,  which  is  very  close  to  the  electron  mobility  of  Pb^j  g^gEug  042^6- 

For  Pbp  527^00  QygTe/PbTe  MQW  samples,  we  found  high  electron  mobilities  at  large  barrier  thick¬ 
ness,  confirming  our  expectations  of  2D  conduction  for  these  wide  barriers.  To  investigate  the  theoretical 
predictions  [2], [10],  various  Pfag  927EUQ o^gTe/PbTe  MQW  samples  were  grown  with  PbTe  well  widths  in 
the  range  15-50  A  and  Pbp  927EUQ  073Te  barrier  widths  in  the  range  360-540  A.  The  Eu  content  in  the 
barriers  was  kept  at  x  =  0.073  +  0.004,  which  gives  a  difference  in  band  gap  between  the  PbTe  wells  and 
the  Pbp  927EUQ  Q23Te  barriers  of  ~  0.4  eV  at  300  K.  The  published  works  agree  that  a  type  I  band 
alignment  occurs  for  a  Pbj_^u^Te/PbTe  superlattice  in  the  range  of  compositions  of  x  smaller  than  0.10. 
Since  for  x  <  0.05  the  band  offsets  are  nearly  symmetric  [11], [12],  we  would  expect  for  x  =  0.073  that 
the  conduction  band  offset  or  barrier  height  is  ~  0.2  eV.  Modulation-doped  Pb^  953EUQ  Q^^Te/PbTe  MQW 
samples  with  very  wide  barrier  layers  and  low  carrier  concentrations  have  been  shown  to  exhibit  en¬ 
hanced  carrier  mobilities  at  77  K  because  of  the  spatial  separation  of  electrons  from  their  donors  [5].  In 
modulation  doping  only  the  wide-gap  barrier  layer  material  is  doped  leaving  the  PbTe  quantum-well  layer 
relatively  free  of  impurities.  Electrons  are  transferred  laterally  from  the  donor  states  of  the  higher-lying 
conduction  band  edge  of  the  Pbg  927EUQ  073Te  barriers  to  the  lower  conduction  band  states  in  the  PbTe 
well  layer.  Thus,  a  spatial  separation  of  electrons  from  the  impurities  occurs.  All  of  our  MQW  samples 
were  modulation  Bi  doped. 

The  measured  thermoelectric  parameters  of  the  Pbg  927EUQ  g73Te/PbTe  MQW  samples,  including 
carrier  concentrations,  carrier  mobilities,  Seebeck  coefficients,  and  thermoelectric  power  factors,  are 
given  in  Table  3-1.  The  Seebeck  coefficients  are  also  displayed  in  Figure  3-2  as  a  function  of  carrier 
concentration  along  with  values  for  our  single-layer  PbTe  films  and  numerous  selected  high-quality  bulk 
n-type  PbTe  samples  [13].  At  low  carrier  concentrations  the  average  300-K  carrier  mobility  for  the  bulk 
samples  was  1600  cm^A^  s,  which  decreased  to  940  cm^A^  s  at  a  carrier  concentration  of  3.9  x 
10*^  cm“^.  The  S  vs  n  data  for  the  bulk  PbTe  was  fit  by  a  straight  line  given  by  the  simple  expression 

S  OuV/K)  =  -477  +  175  logj0(n/lOi7  cm-3)  .  (3.6) 

From  this  expression  and  the  slow  mobility  change  with  carrier  concentration,  we  get  a  maximum 
thermoelectric  power  factor  of  38  /iW/cm  at  ~  6  x  10^®  cm~^,  which  is  consistent  with  the  best 
reported  values  for  PbTe.  Using  k)  =  20  mW/cm  K,  we  calculate  =  0.45. 
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Our  single-layer  PbTe  films  had  somewhat  lower  electron  mobilities  than  high-quality  bulk  PbTe. 
This  is  presumably  a  result  of  unoptimized  growth.  However,  one  of  our  MQW  samples  had  a  very  high 
mobility  (1420  cmW  s) — higher  than  that  of  our  single-layer  PbTe  films  and  at  least  equal  to  the 
mobility  for  the  best  bulk  PbTe  samples.  That  this  is  achieved  when  the  MQW  film  is  composed  mostly 
(95%)  of  very  low  mobility  (45  cm^A^  s  at  300  K)  Pbo.927Euo  o^gTe  material  is  strong  evidence  that  pure 
2D  conduction  is  occurring  and  that  modulation  doping  can  be  effective  in  these  Pbo  927^^0.073^®^^^® 
MQW  samples  with  nearly  optimized  growth  conditions.  There  is  a  wide  variation  in  the  carrier  mobility 
from  the  best  sample  to  the  worst,  and  it  is  believed  that  mobility  is  quite  sensitive  to  the  growth  variables 
such  as  effusion  cell  flux  ratios.  In  general,  our  Pbo  927^^^0.073^®^^^®  films  have  somewhat  lower 

electron  mobilities  than  the  single-layer  PbTe  films,  although  20x  higher  than  that  of  the  Pbo  927EUQ  Q23Te 
barrier  material.  No  focused  effort  was  applied  to  optimization  of  the  epitaxial  growth. 


TABLE  3-1 


Thermoelectric  Properties  of  Pbp  g27EuQ  Q73Te/PbTe  Muitiple-Quantum-Well  Structures 


Sample 

No. 

Thickness 
of  PbTe 
Wells 
(A) 

Carrier 

Concentration 

(cm“®) 

Carrier 
Mobility 
(cm2/V  s) 

Seebeck 

Coefficient 

(AiV/K) 

2D  Power 
Factor 

P/=2D 

(mV/cm  K^) 

Calculated 
Figure  of 
Merit 

^20’’ 

T-225 

19 

1.1 

X 

10^9 

1420 

-225 

130 

1.23 

T-229 

20 

1.1 

X 

10^9 

710 

-235 

66 

0.78 

T-235 

23 

5.9 

X 

10^8 

650 

-303 

55 

0.73 

T-246 

15.5 

7.5 

X 

10^8 

1050 

-222 

62 

0.72 

T-236 

23 

6.7 

X 

1018 

640 

-282 

54 

0.71 

T-218 

45 

4.8 

X 

10I8 

950 

-174 

22 

0.29 

T-222 

34 

4.8 

X 

10I8 

890 

-166 

19 

0.24 

T-208 

40 

3.7 

X 

1018 

800 

-188 

17 

0.23 

T-215 

41 

5.1 

X 

1018 

850 

-157 

17 

0.22 

The  Seebeck  coefficient  values  of  our  PbTe  epilayers  are  equivalent  to  values  of  the  bulk  PbTe  and  nicely 
fill  in  the  data  gap  in  Figure  3-2  near  6  x  10^*  cm'^.  The  Seebeck  coefficient  values  of  the  MQW  films 
show  a  decrease  with  increasing  carrier  concentration  similar  to  that  of  bulk  material.  In  the  theoretical 
calculations  [2],[10]  of  Z^^T,  carrier  mobihty  and  thermal  conductivity  in  the  quantum  wells  were 
assumed  equal  to  bulk  values  and  independent  of  well  thickness.  Thus,  from  Equation  (3.4)  we  need  only 
examine  ^n,  and  since  S  varies  with  n,  we  can  simply  compare  quantum-well  S  values  with  bulk 
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Figure  3-2.  Seebeck  coefficient  vi  carrier  concentration  for  n-type  bulk  PbTe,  epitaxial  PbTe,  and 
Pbg  ^27^^0  073^^^^^'^^  wells. 


values  near  the  optimum  carrier  concentration  (-6x10^^  cm“^).  Note  that  the  samples  with  ~  40- A  well 
widths  have  S  values  very  close  to  the  bulk  curve.  As  the  well  width  was  decreased  S  increased,  and  at 
well  widths  near  20  A,  we  obtained  Seebeck  coefficients  about  2x  higher  than  those  obtained  with  high- 
quality  bulk  material.  This  rapid  increase  with  decreasing  well  width  is  in  agreement  with  the  calculations 
[2], [10].  Well  widths  thinner  than  17  A  did  not  result  in  improvements,  possibly  because  of  interface 
effects  or  because  2D  conduction  is  no  longer  achieved.  The  thickness  of  these  thin  wells  corresponds 
to  only  four  to  six  atom  layers  or  two  to  three  molecular  layers. 

The  thermoelectric  power  factors  {S^a  and  are  presented  in  Figure  3-3  as  a  function  of  the 

carrier  eoncentration  at  300  K.  The  dotted  curve  was  calculated  by  fitting  the  bulk  Seebeck  coefficient 
and  mobility  data,  and  provides  information  near  the  optimum  bulk  carrier  concentration  of  ~  6  x 
10^*  cm“^,  where  ~  38  jUW/cm  is  expected.  Because  of  slightly  lower  electron  mobilities,  the 
power  factors  for  our  single-layer  PbTe  films  are  below  the  bulk  PbTe  curve.  The  2D  power  factors  S^a.^ 
for  thick-quantum-well  films  are  somewhat  further  below  the  bulk  PbTe  curve,  mainly  because  of  low 
mobility.  However,  for  the  thin  PbTe  quantum  wells,  significantly  higher  2D  power  factors  are  seen 
despite  their  generally  lower  mobilities.  For  the  one  high-mobility  sample  T-225  the  power  factor  is 
130  jUW/cm  K^,  the  highest  reported  to  date  for  good  thermoelectric  materials. 
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Figure  3-3.  Thermoelectric  power  factor  vs  carrier  concentration  for  n-type  bulk  PbTe,  epitaxial  PbTe,  and 
PbQg27^^0  073^e/PbTe  quantum  wells. 


The  2D  figures  of  merit  for  the  PbTe  quantum  wells  shown  in  Table  3-1  were  calculated  using 
Equations  (3.3)  and  (3.4).  For  all  five  narrow-quantum-well  samples,  exceeded  the  expected 
=  0.45  of  bulk  PbTe.  For  sample  T-225,  the  combined  enhancement  of  the  Seebeck  coefficient  and 
mobility  resulted  in  Z.^^T  =  1.23,  a  result  that  breaks  the  decades-old  ZT  =  1.0  barrier  for  thermoelectric 
materials  at  300  K.  These  results  establish  the  proof  of  principle  of  an  enhanced  ZT  at  300  K  in  thermo¬ 
electric  quantum  wells,  presumably  because  of  the  large  density  of  states  per  unit  volume  increase  that 
occurs  for  small  well  widths  in  a  quasi-2D  geometry.  Seebeck  coefficients  are  2x  larger  than  those  of 
high-quality  bulk  PbTe  for  quantum-well  widths  in  the  range  17-23  A.  The  combined  enhancements  of 
the  Seebeek  coefficient  and  mobility  resulted  in  an  increase  in  the  300-K  ZT  of  PbTe  from  -  0.45  for 
bulk  material  to  Z^^T  >  1.2  in  the  best  quantum  wells.  Although  much  remains  to  be  done,  as  the  best 
overall  or  3D  ZT  of  these  MQW  structures  was  only  0.11,  this  work  represents  a  first  step  toward  the 
goal  of  achieving  materials  with  3D  ZT  >  1 .0. 

T.  C.  Harman  R.  E.  Fahey 
D.  L.  Spears  M.  P.  Walsh 
M.  J.  Manffa 
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4.  SUBMICROMETER  TECHNOLOGY 


4.1  PLASMA-DEPOSITED  SILYLATION  RESIST  FOR  193-nm  LITHOGRAPHY 

Dry-deposited,  dry-developed  photoresists  hold  numerous  potential  advantages  over  conventionally 
synthesized  spin-applied,  wet-developed  photoresists.  All-dry  resists  enable  lithography  on  materials  that 
may  react  with  solvents  or  liquid  developers,  those  which  require  minimal  handling,  or  those  prone  to 
contamination  or  oxidation.  Photoresists  deposited  by  plasma-enhanced  chemical  vapor  deposition  (PECVD) 
can  be  used  on  any  size  or  shape  substrate,  over  varying  degrees  of  topography,  and  permit  process 
flexibility  through  in-situ  synthesis  tailored  to  a  specific  application.  They  also  avoid  costly  waste  and 
disposal  problems  associated  with  conventional  spin-on  resists.  Our  initial  work  on  negative-tone  PECVD 
organosilicon  thin  films  for  193-nm  lithography  [1]  and  more  recent  work  by  Wiedman  and  Joshi  [2]  on 
organosilicon  imaging  layers  for  mid-ultraviolet  lithography  show  the  potential  viability  of  all-dry 
photoresist. 

Plasma-deposited  photoresist  synthesis  has  different  constraints  than  conventional  polymer  synthe¬ 
sis.  An  advantage  is  that  the  desired  chemical  properties  of  the  polymer  are  not  limited  by  the  precursor 
materials  but  can  be  created  through  plasma  synthesis.  A  drawback  to  this  synthesis  method  is  that  as  a 
single-wafer  process,  reproducibility  is  essential.  Also,  the  deposition  rate  must  be  at  least  on  the  order 
of  0.3  ;Um/min.  In  addition,  plasma  polymerization  has  sharply  unstable  regimes  with  respect  to  deposi¬ 
tion  of  uniform,  particle-free  polymer  films.  Deposition  variables  such  as  precursor  ratios,  flow  rates, 
deposition  pressure,  power,  electrode  spacing  (geometry  of  reactor),  and  temperature  must  be  carefully 
optimized  to  avoid  gas-phase  polymerization  or  highly  nonuniform  deposition. 

In  our  experiments,  all  films  were  deposited  from  liquid  hydrocarbon  precursors  and  oxygen-based 
gases  in  a  conventional  13.56-MHz  RF-powered  parallel-plate  reactor.  Gases  were  introduced  through  a 
water-cooled  shower  head  in  the  top  electrode  at  flow  rates  from  1  to  400  seem.  Films  were  deposited 
at  relatively  high  pressures  (300-3000  mTorr)  and  various  RF  powers  (50-400  W)  onto  4-in.  wafers 
robotically  loaded  onto  the  bottom,  temperature-controlled  (powered)  electrode.  A  variety  of  techniques, 
including  ultraviolet-visible  spectroscopy,  Fourier  transform  infrared  spectroscopy  (FTIR),  thermal-me¬ 
chanical  analysis  (TMA),  and  ellipsometry  were  used  to  characterize  the  photoresist.  Patterned  exposures 
were  performed  using  Lincoln  Laboratory’s  0.35-numerical-aperture  modified  193-nm  GCA  lOx  stepper. 
Samples  were  silylated  using  dimethylsilyldimethylamine  (DMSDMA)  at  10  Torr  for  various  times  and 
temperatures  in  a  computer-controlled  silylation  system.  Pattern  transfer  was  performed  with  oxygen 
plasmas  both  in  a  conventional  reactive  ion  etching  system  at  -200-V  self-bias  and  in  a  high-ion-density 
helicon  reactor. 

Our  benchmark  in  developing  an  all-dry  positive-tone  silylation  resist  for  193-nm  lithography  is 
spin-applied  polyvinylphenol,  also  known  as  polyhydroxystyrene  (PHOST).  We  began  plasma  synthesis 
by  depositing  films  from  high-vapor-pressure  unsaturated  cyclic  hydrocarbon  precursors  and  oxygen- 
based  gases,  since  no  volatile  aromatic  compounds  containing  hydroxyl  groups  are  available  as  precursors 
to  synthesize  a  chemical  analog  to  PHOST  using  PECVD.  Figure  4-1  shows  the  FTIR  spectrum  of  a  film 
deposited  from  benzene  and  oxygen  and  for  comparison  the  spectrum  of  an  equal  thickness  film  of 
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Figure  4-1.  Fourier  transform  infrared  spectroscopy  (FTIR)  absorbance  spectra  of  polyhydroxystyrene  and  a  film 
deposited  from  benzene  and  oxygen.  The  bands  are  identified,  which  are  signatures  of  the  O-H,  and  C—0 

moieties. 


PHOST.  Note  that  both  films  have  a  characteristic  broad  hydroxyl  peak  centered  around  3400  cm"i  and 
both  have  an  aliphatic  C-H^  stretch  centered  around  2950  cm"!  and  aromatic-related  peaks  at  1600  and 
3050  cm'^  The  existence  of  an  aromatic  component  in  the  polymer  films  is  desirable  for  etch  resistance 
in  subsequent  patterning  steps.  Beyond  the  expected  peak  broadening  associated  with  the  amorphous 
nature  of  the  PECVD  film,  the  plasma-deposited  polymer  has  a  carbonyl  peak  at  around  1700  cwT^  not 
found  in  PHOST. 

Among  the  many  variables  in  plasma  synthesis,  we  have  found  that  the  deposition  pressure,  RF 
power,  and  oxygen  to  hydrocarbon  flow  ratio  have  the  most  dramatic  effect  on  deposition  rate,  glass 
transition  temperature  (T^),  and  the  ratio  of  0-H  to  C-H^  in  the  films.  In  general,  the  deposition  rate  is 
most  strongly  affected  by  the  pressure.  All  other  parameters  being  constant,  the  deposition  rate  increases 
almost  linearly  from  150  to  650  nm/min  as  pressure  is  increased  from  500  to  1500  mTorr.  The  of  these 
films,  as  deduced  from  the  onset  of  softening  in  a  TMA  trace,  is  most  strongly  affected  by  the  RF  power. 
For  example,  films  deposited  at  1500  mTorr  at  powers  of  100,  200,  and  300  W  had  a  of  80,  90,  and 
105°C,  respectively. 

A  key  component  of  the  plasma  synthesis  for  this  resist  scheme  is  the  formation  of  hydroxyl  bonds 
in  the  polymer  which,  when  exposed  to  a  silicon-containing  gas,  react  to  incorporate  silicon  into  the  film. 
Figure  4-2  plots  the  deposition  rate  and  the  integrated  FTIR  O-H  to  C-H^  ratio  for  films  deposited  from 
benzene  and  oxygen  for  varying  oxygen  flow  rates.  The  pressure  and  power  were  fixed  at  700  mTorr  and 
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Figure  4-2.  Deposition  rate  and  integrated  [0-H]/[C-HJ  ratio  for  films  deposited  from  benzene  at  various  flow 
rates  of  oxygen. 


150  W,  respectively.  The  benzene  was  introduced  directly,  i.e.,  without  the  use  of  a  mass  flow  controller. 
As  the  oxygen  flow  rate  increases,  the  deposition  rate  decreases  while  the  O-H  to  C-H^  ratio  increases. 
We  believe  that  the  decrease  in  deposition  rate  with  increased  oxygen  flow  is  a  result  of  the  competition 
between  deposition  and  etching  processes.  The  increase  in  O-H  to  C-H.  ratio  with  increased  oxygen  flow 
rate  is  attributed  to  an  increase  in  availability  of  oxygen. 

The  reactivity  of  the  hydroxyl  sites  to  a  silylamine  gas  is  shown  in  Figure  4-3.  This  figure  plots 
the  FTIR  absorbance  spectra  of  a  film  deposited  from  benzene  and  oxygen  after  exposure  to  DMSDMA 
at  55°C  and  10  Torr  for  various  times.  As  the  silylation  time  increases,  more  silicon  is  incorporated,  as 
indicated  by  the  increase  in  the  Si-H  stretch  peak  at  2185  cm“^  and  in  the  C-H3  asymmetric  stretch  peak 
at  2950  cm“^  both  of  which  are  associated  with  a  dimethylsilyl  group.  There  is  also  a  corresponding 
decrease  in  the  O-H  peak  and  an  increase  in  the  Si-O-C  (located  at  923  cm“*  but  not  shown)  with  longer 
silylation  time,  indicative  of  the  reaction  between  the  silicon-containing  gas  and  the  hydroxyl  groups  [3]. 
The  silylation  process  for  these  films  is  linear  as  seen  both  with  FTIR  and  scanning  electron  microscopy 
cross-sectional  staining,  indicating  that  the  rate-limiting  step  is  not  Fickian  diffusion.  The  absolute  silylation 
rate  is  a  strong  function  of  temperature;  the  films  exhibit  activation  energies  for  diffusion  from  0.55  eV, 
which  is  similar  to  that  reported  for  DMSDMA  in  PHOST  [4],  to  >  1.0  eV,  depending  upon  deposition 
conditions. 
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Figure  4-3.  Time  sequence  ofFTIR  absorbance  spectra  of  a  film  deposited  from  benzene  and  oxygen  silylated  with 
dimethylsilyldimethylamine  (DMSDMA)  at  45  °C  and  10  Torrfor  silylation  times  of  30,  45,  60,  and  120  s. 


The  photosensitivity  of  these  PECVD  films  is  a  direct  measure  of  how  effective  the  material  is, 
when  crosslinked  by  6.4-eV  photons,  in  retarding  the  diffusion  of  a  gaseous  silylating  agent.  For  many 
spin-applied  phenolic  resists,  193-nm  radiation  is  strongly  absorbed  causing  homolytic,  free-radical 
crosslinking  of  the  resist  [4].  The  presence  of  a  cyclic/aromatic  component  in  these  PECVD  films  (see 
Figure  4-1)  is  believed  to  be  important  in  facilitating  the  crosslinking  step,  thereby  controlling  photosen¬ 
sitivity.  We  have  used  the  sequence  of  open  frame  exposure,  silylation  to  a  fixed  depth,  and  pattern 
transfer  with  10%  overetch  to  assess  relative  photospeeds  of  these  resists.  We  have  seen  photospeeds 
from  30  to  300  mJ/cm^  for  polymers  deposited  from  different  precursors  and  under  different  plasma 
conditions,  compared  to  45  mJ/cm^  for  PHOST  silylated  to  a  thickness  of  200  nm. 

The  lithographic  contrast  of  the  silylation  process  is  evaluated  by  inferring  the  silicon  uptake  from 
the  integrated  Si-H  peak  as  a  function  of  dose.  As  in  PHOST,  the  silylation  contrast  is  rather  low,  ~  1, 
but  the  contrast  after  pattern  transfer  is  much  higher,  approaching  20.  The  optical  density  (OD)  at  193 
nm  is  also  important  to  the  lithographic  performance  of  this  resist,  because  it  determines  the 
photocrosslinking  depth.  For  spin-applied  193-nm  silylation  resists,  the  crosslinked  layer  is  typically 
<  70  nm  thick,  while  the  depth  of  silicon  incorporation  necessary  for  pattern  transfer  is  much  greater.  This 
relatively  thin  crosslinked  layer  does  not  prevent  isotropic  diffusion  under  it,  thus  potentially  compromis¬ 
ing  the  ultimate  lithographic  performance  for  sub-0.25-;Um  lithography.  The  OD  for  films  formed  from 
benzene  and  oxygen  is  2-4  /mr^,  much  less  than  that  of  PHOST  (25  and  corresponding  to  the 

depth  necessary  to  prevent  isotropic  diffusion.  Selective  silylation,  and  thus  the  lithographic  capability, 
of  PECVD  silylation  photoresist  is  demonstrated  in  Figure  4-4,  which  shows  an  example  of  residue-free 
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Figure  4-4.  Scanning  electron  micrograph  of  0.4-fm  features  printed  with  the  positive-tone  all-dry  silylation  resist 
deposited  from  cycloheptatriene  and  oxygen.  Exposure  was  done  on  a  modified  193-nm  GCA  lOx  stepper  at  a  dose 
of  135  mJ/cm^.  Silylation  was  performed  using  DMSDMA  at  85  °C  and  10  Torr  for  90  s.  Pattern  transfer  was 
performed  in  a  helicon  etcher  using  oxygen  for  80  s  yielding  200%  overetch. 


0.4-fJm  features  patterned  in  a  0.70-jUm- thick  film  deposited  from  cycloheptatriene  and  oxygen.  The 
sample  was  exposed  to  a  dose  of  135  mJ/cm^,  silylated  at  85°C  and  10  Torr  with  DMSDMA  for  90  s, 
and  pattern  transferred  for  80  s  in  a  helicon  reactor  yielding  200%  overetch  (oxygen  plasma).  These  are 
initial  results  of  this  novel  process.  We  expect  that  significantly  higher  resolution  can  be  obtained  as  the 
process  is  further  refined,  and  with  the  use  of  higher-numerical-aperture  projection  optics. 
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FIELD  EMISSION  FROM  NITROGEN-DOPED  DIAMOND 

Diamond  has  a  negative  electron  affinity  (NEA)  [5]— [8],  that  is,  the  minimum  energy  of  electrons 
in  the  conduction  band  is  above  the  minimum  energy  of  electrons  in  vacuum.  This  has  made  possible  the 
realization  of  electron  diode  emitters  [9],[10],  and  has  motivated  others  to  characterize  electric  field 
emission  from  diamond  [11]-[14].  However,  to  take  advantage  of  diamond’s  NEA  property,  a  source  of 
electrons  must  be  present  either  in  the  conduction  band  or  in  a  sub-band  near  the  vacuum  energy  level 
[14]. 


Nitrogen  forms  a  deep  electron  donor  ~  1.7  eV  below  diamond’s  conduction  band  [15],  and  we  have 
observed  electron  emission  from  nitrogen-doped  diamond  at  vacuum  electric  fields  of  0-1-V  /rni"'. 
However,  because  these  diamonds  are  very  resistive,  a  substantial  voltage,  several  thousand  volts,  is 
dropped  in  the  1-  to  2-mm-thick  diamond  substrate.  These  results  should  be  contrasted  to  the  behavior 
of  polycrystalline  boron-doped  (p-type)  diamond.  This  material  is  highly  conductive  and  does  not  exhibit 
a  significant  voltage  drop  across  the  substrate,  but  it  requires  high  fields,  20-50  V  ^  to  obtain 
emission. 

Type  Ib  diamond  substrates,  doped  to  ~  10^^  cm“^  with  nitrogen,  were  obtained  from  commercial 
sources  for  use  in  this  study.  One  set  of  experiments  used  (lOO)-oriented,  cut  and  polished  diamonds.  The 
other  experiments  used  wire-drawing  dies  (synthetic  high-pressure  diamonds  used  as  dies  for  wire  draw¬ 
ing)  that  consisted  of  cubo-octahedrons,  faceted  on  (111)  and  (100)  crystal  planes,  2-3  mm  on  a  side.  All 
these  substrates  were  cleaned  [14]  prior  to  mounting  in  a  10"^  Pa  ion-  or  a  cryo-pumped  vacuum  system. 

Electrical  contact  to  the  boron-doped  diamond  can  be  made  with  a  clip,  but  type  Ib  substrates  are 
insulating  and  the  back-metal  contact  was  found  to  be  the  limiting  factor  for  both  conduction  through  and 
emission  from  the  diamond.  Contact  to  the  diamond  was  made  by  enhancing  field  emission  of  electrons 
into  the  diamond  with  a  rough  metal-diamond  interface.  Two  techniques  were  employed  to  obtain  this 
interface.  One  used  either  Ni  or  Ni  salts  on  the  back  of  the  diamond,  annealed  in  a  mixture  of  Ar  and 
H  at  900-1080°C  for  several  minutes.  The  Ni  or  the  Ni  metal  formed  from  the  metal  salts  during 
annealing  etches  pits  in  the  substrate  by  catalyzing  a  reaction  between  diamond  and  H2  to  form  CH4.  The 
other  technique  used  anisotropic  dry  etching  [16]  to  roughen  the  diamond  surface,  which  was  later 
covered  with  vacuum-evaporated  Ni  to  make  electrical  contact.  Both  techniques  increased  the  conductiv¬ 
ity  though  the  diamond.  Figure  4-5(a)  compares  the  conductivity  through  1  mm  of  type  lb  diamond  for 
a  smooth  and  rough  interface.  Conduction  through  these  diamonds  usually  approximates  a  straight  line 
when  plotted  in  Fowler-Nordheim  [17]  coordinates  (current  through  the  diamond  divided  by  the  voltage 
across  the  diamond  squared  plotted  against  the  inverse  of  the  voltage),  as  shown  by  Figure  4-5(b). 

Field  emission  measurements  were  made  by  placing  a  l-nun-diam  Mo  anode  at  various  heights 
above  the  sample.  As  seen  in  Figure  4-6,  the  emission  is  nearly  independent  of  substrate-anode  spacing 
from  touching  to  400  jum,  and  only  decreases  when  this  spacing  is  larger  than  the  lateral  dimensions  of 
the  diamond  substrate,  ~  5  mm.  An  estimate  of  the  field  between  the  substrate  and  the  anode  can  be 
obtained  by  plotting  the  voltage  required  to  obtain  a  given  emission  current  for  several  substrate-anode 
spacings.  The  slope  of  the  straight  line  fit  to  the  data  is  the  average  vacuum  field,  and  the  intercept  at 
zero  spacing  is  the  internal  potential  across  the  substrate.  Figure  4-7  compares  data  from  nitrogen-doped 
diamond  and  boron-doped  diamond  when  plotted  in  this  manner.  Conductive  boron-doped  diamond  does 
not  support  a  measurable  potential  across  the  substrate,  but  requires  a  substantial  field  to  obtain  emission. 
Type  Ib  diamond  requires  a  vacuum  field  of  only  ~  0.2  V  /lm-^  but  a  potential  of  several  kilovolts 
appears  across  the  substrate. 
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VOLTAGE  ACROSS  DIAMOND  (V) 

(a) 


(b) 


Figure  4-5.  (a)  Comparison  of  conductivity  through  two  halves  of  the  same  1-mm-thick,  polished,  type  Ib  diamond. 
The  back  contact  consisted  of  a  smooth  interface  with  100  nm  of  electron-beam-evaporated  Ni  and  a  roughened 
interface  with  a  Ni  salt  deposited  on  the  diamond  and  annealed  in  H2  and  Ar  at  1080°C  for  30  min  (an  additional 
100  nm  ofNi  was  evaporated  on  the  back  of  the  diamond  to  ensure  electrical  contact),  (b)  Results  for  the  roughened 
interface  plotted  in  Fowler-Nordheim  coordinates. 
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ANODE  VOLTAGE  (V) 
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Figure  4-6.  Emitted  current  as  a  function  of  anode  voltage  for  several  substrate-anode  spacings.  The  back  contact 
was  roughened  as  described  in  the  caption  of  Figure  4-5. 
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Figure  4-7.  Voltage  required  for  field  emission  of  1  X  I0~^  A  cm  ^  for  boron-doped  diamond  {open  circles)  and 
type  Ib  diamond  (solid  circles)  as  a  function  of  substrate-anode  spacing. 


The  emission  from  both  polished  and  as-grown  wire-drawing  type  Ib  diamonds  appears  the  same. 
The  electrons  are  either  emitted  as  a  diffuse  unstructured  beam  or  as  many  small  beamlets.  Diffuse  beams 
are  emitted  from  diamonds  that  have  been  coated  with  Cs,  as  described  in  [14],  while  the  uncoated 
diamonds  usually  emit  beamlets.  These  beamlets  are  emitted  at  angles  up  to  45°  from  the  substrate 
normal.  By  using  a  movable  phosphor  screen  it  was  determined  that  small  regions  <  100  fim  in  diameter 
were  responsible  for  several  beamlets  emerging  at  angles  several  tens  of  degrees  from  each  other.  When 
the  emitted  current  is  >  0.1  mA,  glowing  channels  often  appear  in  the  diamond  that  extend  from  the  back 
contact  to  the  top  of  the  diamond.  The  source  of  beamlets  appears  to  correspond  to  the  top  of  these 
glowing  channels.  Often  the  beamlets  would  be  so  energetic  that  they  would  miss  the  anode  ~  100  jum 
from  the  diamond  surface  and  strike  the  anode’s  metal  fixturing  several  millimeters  away.  When  the 
emitted  current  is  ~  0.1-1  mA,  the  electrons  emitted  from  the  diamond  have  sufficient  energy  to  cause 
the  anode  to  glow  red  hot.  This  is  known  to  require  power  levels  >  1  W,  which  necessitates  that  the 
average  electron  energy  be  several  kilovolts.  This  electron  energy  is  considerably  larger  than  the 
~  100  eV  that  the  electrons  could  obtain  traveling  from  the  diamond  surface  to  the  anode  when  the 
vacuum  field  is  <  1  V  /tm~*. 

The  emission  characteristics  of  the  boron-doped  and  nitrogen-doped  diamond  can  be  explained  as 
follows.  Boron-doped  diamond  has  effectively  no  conduction  band  electrons  and  a  work  function  of 
4.8-7. 8  eV  [7],  requiring  electric  fields  of  20-50  V  firrr^  to  obtain  electron  emission  from  rough  poly¬ 
crystalline  surfaces,  the  same  field  necessary  for  field  emission  from  high-work-function  metals,  like  W 
or  Mo.  For  nitrogen-doped  diamond,  where  the  substrate  is  insulating  at  low  fields,  <  0.1  V  jum“^  the 
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emission  is  dependent  upon  the  back  interface.  The  nitrogen  dopant  forms  a  depletion  region  at  the  back 
metal-diamond  interface  which  allows  for  tunneling  of  electrons  into  the  diamond.  A  rough  surface  at  the 
interface  enhances  the  electric  field  locally  and  facilitates  electron  tunneling.  Once  in  the  diamond, 
electrons  travel  through  the  diamond  under  the  influence  of  an  internal  electric  field  and  are  emitted  from 
the  diamond  surface.  The  NEA  property  allows  the  emission  of  electrons  into  the  vacuum  with  minimal 
vacuum  field.  Some  of  the  emitted  electrons  appear  to  leave  the  diamond  surface  with  considerable 
energy,  >  1  keV.  We  are  uncertain  where  the  electrons  are  accelerated.  This  energy  could  be  due  to  high 
fields  just  above  or  just  below  the  diamond  surface  or  to  electron  channeling  through  the  diamond.  The 
diamonds  often  glow  during  emission,  and  this  may  indicate  the  presence  of  conduction  channels  [18] 
which  can  form  in  semiconductors  under  high  internal  fields.  Where  these  several  kilovolt  electrons  are 
accelerated  is  being  investigated.  Several  groups  have  reported  electrons  being  accelerated  to  high 
energies,  5-60  eV,  in  insulators  [19]  and  semiconductors  [20].  Guttler  [21]  has  predicted  higher  electron 
energies,  ~  1  keV,  in  diamond  when  subjected  to  fields  >  100  V  /tm~^ 
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5.  mGH  SPEED  ELECTRONICS 


5.1  SELF-ALIGNED  COMPLEMENTARY  GaAs  MISFET  USING  A 
LOW-TEMPERATURE-GROWN  GaAs  GATE  INSULATOR 

A  self-aligned  «-channel  metal-insulator-semiconductor  FET  (MISFET)  with  very  low  gate  leakage 
current  has  been  developed  using  a  low-temperature-grown  (LTG)  GaAs  gate  insulator  [1].  The  MISFET 
is  ideal  for  digital  circuits  because  it  provides  high  noise  margin  and  low  static  power  consumption. 
Recently  developed  self-aligned  p-  and  n-channel  MISFETs  fabricated  on  the  same  wafer  are  reported 
here  to  demonstrate  the  feasibility  of  using  LTG-GaAs  MISFETs  for  complementary  logic  circuits. 

The  epitaxial  layers  were  grown  by  molecular  beam  epitaxy  on  a  semi-insulating  GaAs  substrate, 
as  shown  in  Figure  5-l(a).  The  layers  consist  of  a  3000-A-thick  GaAs  channel,  a  100-A-thick  AlAs 
barrier,  a  300-A-thick  LTG-GaAs  gate  insulator,  a  100-A-thick  AlAs  barrier,  and  a  100-A-thick  GaAs 
cap.  None  of  the  layers  was  intentionally  doped,  but  the  undoped  GaAs  channel  was  slightly  p-type  with 
a  background  hole  concentration  of  ~  10^^  cm“^.  The  growth  temperature  was  600°C  for  all  the  layers 
except  for  the  LTG  GaAs,  which  was  grown  at  200°C  and  then  annealed  at  600°C  for  10  min  before 
growing  the  AlAs  barrier  and  GaAs  cap  on  top.  The  3-in.  wafer  with  these  epitaxial  layers  was  then 
cleaved  into  several  pieces  for  the  fabrication  of  normally-off  p-channel  and  both  normally-off  and 
normally-on  n-channel  MISFETs.  Each  piece  had  only  one  type  of  MISFET  on  it  and  the  same  mask  set 
was  used  for  all  devices. 

A  simplified  illustration  of  the  fabrication  sequence  is  shown  in  Figure  5- 1(b).  For  the  normally- 
on  n-channel  MISFET,  the  doped  n-channel  was  first  formed  by  Si  implantation  without  a  mask.  The 
implant  energy  and  dose  were  115  keV  and  7  x  10^^  cm~^,  respectively.  Then  the  gate  metal,  consisting 
of  700  A  of  W  and  2000  A  of  Au,  was  deposited  by  evaporation  and  patterned  by  optical  lithography 
and  liftoff  on  all  three  pieces.  Next,  the  source  and  the  drain  regions  were  implanted  using  the  W/Au  gate 
metallization  as  a  self-aligned  mask.  For  the  two  n-channel  pieces,  Si  was  used  to  form  the  regions 
with  an  implant  energy  of  60  keV  and  a  dose  of  3.0  x  10^^  cm“^.  For  the  p-channel  piece,  Be  was  used 
to  form  the  regions  with  an  implant  energy  of  24  keV  and  a  dose  of  4.7  x  10^^  cm“^.  After  the 
implantation,  the  Au  metal  was  etched  off  and  the  wafers  were  annealed  at  825°C  for  10  s  using  a  GaAs 
substrate  as  the  proximity  anneal  cap.  Next,  the  source  and  drain  ohmic  regions  were  defined  in  photo¬ 
resist  and  the  gate  insulator  layers  were  wet  etched  in  the  exposed  ohmic  regions  to  the  top  of  the  GaAs 
channel.  Then  the  ohmic-contact  metal  was  evaporated  and  defined  by  liftoff  using  the  same  photoresist 
that  was  patterned  for  the  etching,  and  this  was  followed  by  alloying.  The  ohmic  metallizations  were  Ni/ 
Ge/Au  for  the  n-channel  and  Au/Zn/Au  for  the  p-channel  MISFETs.  Finally,  devices  were  isolated  by 
proton  implantation.  The  gate  length  and  width  of  the  MISFET  were  1.5  and  100  pm,  respectively,  and 
the  spacing  between  the  source  and  drain  ohmic  contacts  was  7  pm. 

The  drain  current  vs  voltage  characteristics  of  MISFETs  are  shown  in  Figure  5-2.  As 

expected,  the  maximum  of  the  normally-off  n-channel  MISFET  is  small  (40  mA/mm  at  the  gate 
voltage  Vgg  =  +4.5  V),  and  the  device  has  a  very  high  output  resistance.  As  shown  in  Figure  5-2(b),  the 
threshold  voltage  is  1.3  V.  For  the  normally-on  n-channel  MISFET,  the  maximum  is  ~  500  mA/mm 
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because  of  the  Si  implant  in  the  channel.  As  shown  in  Figure  5-2(a),  a  of  —2.4  V  is  required  to  pinch 
off  the  device.  Figure  5-2(c)  shows  that  for  the  p-channel  MISFET  reaches  a  maximum  of  120  mAJ 
mm  at  F  =  -8  V.  The  maximum  of  22  mS/mm  measured  at  =  -5  V  is  comparable  to  that  of  the 
best  GaAs/AlGaAs  p-channel  heterojunction  FET  having  a  similar  gate  length.  The  important  device 
parameters  are  summarized  in  Table  5-1. 
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100  AGaAs 
iooAaias 
300  A  LTG  GaAs 
iooAaias 
3000  A  GaAs 


(1) 


(2) 


NORMALLY-ON  NORMALLY-OFF  NORMALLY-OFF 
n-CHANNEL  n-CHANNEL  p-CHANNEL 

(b) 


Figure  5-1.  (a)  Layer  structure  and  (b)  fabrication  sequence  for  p-  and  n-channel  low-temperature-grown  (LTG) 
GaAs  metal-insulator-semiconductor  FETs  (MISFETs). 
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TABLE  5-1 


Device  Parameters  for  Metal-lnsulator-Semiconductor  FETs 


n-Channel 

Normally-On 

n-Channel 

Normally-Off 

p-Channel 

Normally-Off 

Maximum  /^g  (mA/mm) 

500 

40 

120 

1/gg  for  maximum  1^^  (V) 

+4.0 

+4.5 

-8.0 

Maximum  (mS/mm) 

114 

25 

22 

Threshold  voltage  Vf  (V) 

-2.4 

+1.3 

-1.2 

Reverse  gate  breakdown  (V) 

-10 

^0 

+7 

The  amount  of  the  forward-biased  gate  current  of  all  three  MISFETs  is  very  similar  and  the 
I  -V  relationship  follows  the  typical  trap-limited  current  patterns,  indicating  that  the  conduction  is 
determined  mainly  by  the  LTG-GaAs  layer.  At  1  V  of  forward  bias,  the  gate  current  is  ~  0.1  |iA/mm  gate 
width,  which  is  equivalent  to  6.6  x  10“^*^  A/jUm^  gate  area.  By  using  the  conventional  definition  of 
at  1  mA/mm,  the  gate  forward  tum-on  voltage  is  close  to  3  V,  which  is  much  higher  than  the  0.7  V  for 
a  MESFET.  Since  the  forward-biased  gate  current  is  very  sensitive  to  the  quality  of  the  LTG  insulator, 
our  results  indicate  that  implanting  Si  through  the  LTG-GaAs  gate  insulator  of  the  normally-on  n-channel 
MISFET  did  not  degrade  the  resistivity  of  the  LTG  GaAs. 

The  significance  of  this  work  is  the  demonstration  of  a  new  family  of  MISFETs  ideal  for  comple¬ 
mentary  logic  circuits,  even  though  the  device  parameters  are  not  optimized.  The  results  reported  here 
show  that  the  LTG-GaAs  MISFET  can  sustain  the  ion  implantation  for  both  channel  and  source-drain 
doping  without  noticeable  degradation.  In  addition,  the  threshold  voltage  and  the  drain  current  of  MISFETs 
on  the  same  wafer  can  be  easily  adjusted  by  the  channel  implant.  Other  than  replacing  a  GaAs  substrate 
with  an  epitaxial  wafer,  all  the  processes  for  complementary  LTG-GaAs  MISFETs  are  industry  standards 
for  self-aligned  GaAs  MESFET  integrated  circuits. 

C.  L.  Chen  L.  J.  Mahoney 
K.  B.  Nichols  E.  R.  Brown 
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6.  MICROELECTRONICS 


6.1  NOISE  ANALYSIS  OF  AN  ELECTRONICALLY  SHUTTERED 
AND  BACK-ILLUMINATED  CCD  PIXEL 

Noise  measurements  are  reported  for  an  electronically  shuttered,  back-illuminated  pixel  [1]  struc¬ 
ture  under  high  photon  flux  (shot  noise  limited)  conditions.  Measurements  show  that  over  a  certain 
operating  range,  the  electronic  shutter  suppresses  the  input  shot  noise  more  than  the  corresponding 
average  signal  of  the  photogenerated  electrons,  resulting  in  an  improvement  in  the  overall  signal-to-noise 
ratio.  The  electronic  shutter  causes  noise  suppression  by  correlating  the  collected  photoelectrons  that  are 
generated  by  an  incoming  light  signal  over  an  integration  period.  The  shutter  creates  a  nonlinear  signal- 
dependent  collection  efficiency  similar  to  the  avalanche  photodiode  analyzed  in  [2].  During  an  integration 
period,  this  nonlinear  collection  efficiency  causes  an  anticorrelation  between  the  photoelectrons  already 
accumulated  in  the  pixel  well  and  the  photoelectrons  still  to  be  collected. 

Figure  6-1  shows  impurity  and  potential  profiles  vs  depth  below  the  surface  of  the  collection  well 
where  the  potential  profiles  are  given  for  two  different  bias  conditions.  A  detailed  description  of  the 
electronic  shutter  fabrication  and  operation  are  given  in  [1].  The  bias  conditions  in  Figure  6- 1(c)  cause 
the  depletion  region  to  extend  deep  into  the  high-resistivity  (>  3000  Q  cm)  p-type  substrate.  Photoelec¬ 
trons  created  in  the  substrate  are  accelerated  by  the  electric  field  to  the  pixel  well  as  in  a  conventional 
back-illuminated  charge-coupled-device  (CCD)  pixel.  For  a  constant  photon  flux,  the  collected  electron 
signal  varies  linearly  with  time.  In  Figure  6- 1(b),  the  voltage  bias  on  the  imaging  array  gate  electrode 
causes  the  depletion  edge  to  terminate  in  the  relatively  high  doped  (~  3  x  10^^cm~^)  p  barrier  region.  The 
undepleted  portion  of  the  p  barrier  region  creates  a  potential  barrier  that  reduces  the  number  of  photo¬ 
electrons,  generated  in  the  low-doped  p  substrate,  that  reach  the  pixel  collection  well.  As  more  electrons 
are  collected  in  the  well,  the  depletion  region  continues  to  collapse  and  the  barrier  height  increases.  The 
increasing  barrier  height  with  collected  electron  signal  gives  a  nonlinear  signal-dependent  collection 
efficiency.  As  the  well  continues  to  fill,  the  barrier  height  eventually  becomes  constant  to  the  photoelec¬ 
trons  generated  in  the  bulk. 

Figure  6-2(a)  shows  the  standard  deviation  vs  average  signal  on  a  logarithmic  scale,  and  Fig¬ 
ure  6-2(b)  is  a  plot  of  the  variance  vs  average  signal  on  a  linear  scale.  Each  point  displayed  represents 
one  pixel  in  a  1 6  x  64  section  of  the  64  x  64-pixel  imaging  array.  The  result  for  a  particular  pixel  is  an 
average  calculated  from  400  image  frames.  The  varying  pixel  signals  are  created  by  a  spatially  nonuni¬ 
form  light  source  illuminating  the  array.  At  low  signal  levels,  the  slope  of  the  log  curve  in  Figure  6-2(a) 
is  1/2,  as  would  be  expected  for  linear  collection  of  a  Poisson-distributed  photoelectron  signal.  However, 
as  the  signal  increases,  the  slope  falls  below  1/2,  indicating  noise  suppression  at  the  higher  signal  levels. 
Also  shown  in  Figure  6-2(a)  is  the  pixel  response  that  occurs  when  a  gate  voltage  is  applied  that  causes 
the  depletion  region  to  extend  deep  into  the  substrate.  In  this  case,  the  collection  of  the  photoelectron 
signal  remains  linear  with  time,  as  in  an  unshuttered  pixel,  and  the  slope  of  the  curve  has  a  value  of 
1/2  even  at  the  highest  signal  shown.  Figure  6-2(b)  shows  the  same  tendency  as  in  Figure  6-2(a),  where 
for  a  conventional  shot-noise-limited  signal  the  curve  is  linear,  but  for  the  shuttered  pixel  the  curve  is 
substantially  sublinear.  In  fact,  the  variance  is  almost  a  factor  of  2  lower  than  that  of  the  shot-noise- 
limited  input  signal. 
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Figure  6-1.  (a)  Doping  profile  vs  depth  for  the  shuttered  pixel  collection  well;  (b),(c)  potential  profiles  for  depletion 
region  extending  deep  into  substrate  and  for  depletion  region  terminating  in  highly  doped  p  barrier  region,  respec¬ 
tively. 
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Figure  6-2.  (a)  Standard  deviation  an^  V5  average  number  of  photoelectrons  collected,  shown  on  a  loga¬ 

rithmic  scale;  (b)  same  data  only  now  plotted  as  variance  on^  vs  average  signal  on  a  linear  scale. 
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The  measurements  indicate  that  the  spread  in  the  measured  photoelectron  distribution,  created  by 
a  shot-noise-limited  and  constant  flux  light  signal,  is  less  than  that  corresponding  to  a  Poisson  distribu¬ 
tion.  For  the  highest  signal  level  given  in  Figure  6-2,  the  squared  output  signal-to-noise  ratio  is  increased 
by  a  factor  of  1.75  over  that  of  the  shot-noise-limited  input. 

R.  K.  Reich 
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7.  ANALOG  DEVICE  TECHNOLOGY 


7.1  MICROWAVE  INTERMODULATION  PRODUCTS  IN  YBa2Cu307_^ 

JOSEPHSON  JUNCTIONS  AND  THE  COUPLED-GRAIN  MODEL 

The  generation  of  intermodulation  products  (IMPs)  by  microwave  components  is  undesirable  in 
almost  all  applications.  Normal  metals  and  most  dielectrics  have  linear  electrical  properties  up  to  very 
high  fields  and  currents,  and  therefore  passive  microwave  components  made  from  them  produce  very 
little  intermodulation  distortion.  On  the  other  hand,  components  made  from  the  high-transition-tempera- 
ture  superconductors  have  shown  relatively  higher  nonlinearities,  which  lead  to  power-dependent  losses 
and  higher  levels  of  IMPs  [1].  It  is  believed  that  the  dominant  mechanism  for  these  nonlinearities  is 
Josephson-junction-like  weak  links  between  grains  of  otherwise  nearly  ideal  superconductor  [2].  This 
coupled-grain  model  has  been  used  to  explain  power-dependent  losses  in  superconducting  microwave 
resonators  [3].  The  purpose  of  this  report  is  to  present  a  first  step  toward  determining  whether  the  same 
coupled-grain  model  can  also  explain  the  observed  IMP  generation. 

The  coupled-grain  model  assumes  that  all  the  grains  are  of  equal  size  and  that  all  the  Josephson 
junctions  between  them  are  identical.  In  practice,  a  film  almost  certainly  contains  a  distribution  of  grain 
sizes  and  junction  properties,  which  makes  interpretation  of  experimental  data  difficult.  To  avoid  this 
problem  we  have  tested  our  models  on  an  intentionally  produced  single,  although  spatially  extended, 
Josephson  junction.  These  results  can  easily  be  extended  to  more  complex  coupled-grain  models. 

We  will  analyze  a  single  half-wavelength  resonator  with  a  Josephson  junction  intentionally  placed 
in  the  center  at  z  =  0,  as  shown  in  Figure  7-1.  The  current  density  through  the  junction  is  not  uniform, 
peaking  strongly  at  the  edges  of  the  resonator  [4],  but  each  small  region  of  the  junction  can  be  considered 
as  a  resistively  shunted  Josephson  junction  (RSJ)  with  uniform  current. 


Figure  7-1.  Microwave  resonator  with  Josephson  junction. 
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Currents  applied  simultaneously  through  a  Josephson  junction  at  two  frequencies  and  /2  will 
produce  currents  and  voltages  at  many  orders  of  sum  and  difference  frequencies  because  of  the  junction 
nonlineaiities,  but  the  strongest  IMPs  near  /j  and  /j  will  be  the  third-order  ones  (IM3s),  at  frequencies 
2/J-/2  and  2/2-/1.  We  have  calculated  the  IM3  voltage  produced  across  an  RSJ  by  numerically  solving 
the  usual  RSJ  circuit  equations  [5]  in  the  time  domain.  Figure  7-2  shows  the  normalized  third-order 
intermodulation  voltage  appearing  across  the  junction  vs  the  normalized  current  amplitude.  The  third- 
order  intermodulation  voltage  is  normalized  by  where  is  the  junction  critical  current  and  R  is 
the  junction  shunt  resistance.  The  sum  of  the  two  equal  applied  current  amplitudes  +  I2  =  21  is 
normalized  by  7^^.  The  different  curves  correspond  to  different  values  of  the  reduced  frequency  Q,  given 
by 


Q  = 


2elcR/tt  ’ 


(7.1) 


where  ft)  =  ^(/l+/2)  is  the  average  angular  frequency,  e  is  the  electronic  charge,  and  h  is  Planck’s 
constant  divided  by  2;r. 

Figure  7-2  has  several  interesting  features.  At  low  applied  current  Vj  increases  as  as  would 

be  expected  from  a  third-order  process.  The  amplitude  of  in  this  cubic  region  is  well  approximated 
by  the  expression 


Vi  0.02Q 

f2/f 

f 

/c7?^l  +  2.67Q2+fi4 

UcJ  ' 

UcJ 

(7.2a) 
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Figure  7-2.  Third-order  intermodulation  voltage,  normalized  by  I^R,  vi  the  current  at  each  fundamental  frequency. 
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As  2I/I(^  exceeds  1  for  lightly  damped  (Q  <  1)  junctions,  V3  increases  more  rapidly  than  cubic  and  then 
remains  nearly  constant.  For  Q  «  1,  as  will  be  the  case  for  junctions  studied  here,  we  will  approximate 
this  by  a  step  change  to 


V3 

icR 


=  0.2  , 


>1 . 


(7.2b) 


Figure  7-1  shows  the  center  conductor  of  a  stripline  microwave  resonator  with  a  single  Josephson 
junction  across  the  center.  The  length  of  the  resonator  in  the  propagation  direction  is  an  integral  number 
of  half  wavelengths,  mXj2,  where  is  the  wavelength  at  resonance  on  the  transmission  line  from  which 
the  resonator  is  made,  and  the  mode  number  m  is  a  positive  integer.  The  center  conductor  has  width  w 
in  the  x-direction  and  thickness  b  in  the  y-direction. 

To  find  the  resonator  currents  caused  by  voltages  across  the  Josephson  junction,  we  consider  a  half- 
wavelength-long  resonator  {m  =  1)  with  a  voltage  source  V=V^  coscot  in  the  center,  at  z  =  0,  extending 
across  the  resonator  but  with  negligible  extent  in  the  z-direction.  This  voltage  produces  a  current  in  the 
resonator  dl(x)  =  dl^{x)  cosoot,  where 

dh  (x)  =  bJi  {x)dx  =  VidY{x)  .  (7.3) 


Here  J^{x)  is  the  amplitude  of  the  current  density,  assumed  constant  over  the  resonator  thickness  b.  This 
expression  defines  the  differential  admittance 


WHj)* 

Vi 


(7.4) 


By  reciprocity,  the  short-circuit  current  amplitude  produced  at  z  =  0  by  a  voltage  source  of  ampli¬ 
tude  at  position  x  is 


dl3{x)  =  V3{x)dY{x)  , 


(7.5) 


so  the  total  current  amplitude  from  all  voltage  sources  is 

For  our  purposes  is  the  IM3  voltage  given  by  Equation  (7.2),  so  its  dependence  on  I  has  been  included 
explicitly  in  Equation  (7.6).  Equation  (7.6)  is  the  expression  that  must  be  evaluated  to  calculate  IM3 
currents.  They  can  then  be  compared  to  the  currents  at  the  applied  firequencies.  We  will  first  derive  an 
expression  for  /j(x)/Vj,  then  generalize  Equation  (7.2)  to  yield  V.^{I,x). 

We  seek  a  relation  between  the  resonator  current  and  voltage  amplitudes  when  the  voltage  appears 
at  z  =  0  between  two  halves  of  a  resonator.  This  resonator  is  capacitively  coupled  at  both  ends  to 
transmission  lines  of  characteristic  impedance  Zq.  This  admittance  Fj  can  be  calculated  by  standard 
techniques,  bearing  in  mind  that  at  resonance  the  length  of  each  section  of  transmission  line  differs 
slightly  from  1^/4  to  account  for  the  capacitive  coupling  to  Zq.  The  result  at  resonance  is  approximately 
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(7.7) 


Yi= - f - - - ^  , 

2Zo\aXg  +  {a)CZoY 

where  C  is  the  capacitance  at  each  end  of  the  resonator  coupling  it  to  external  transmission  lines  termi¬ 
nated  in  their  characteristic  impedance  Zq,  ®is  the  angular  frequency,  and  a  is  the  attenuation  coefficient 
of  the  resonator  transmission  line.  The  characteristic  impedance  of  the  resonator  is  also  Zq.  The  derivation 
of  Equation  (7.7)  has  assumed  aX^,  coCZq  «  1. 

Van  Duzer  and  Turner  [6]  give  an  approximate  analytic  form  for  the  eurrent  density  distribution 
of  an  isolated  superconducting  strip  when  wb»X^,  which  is  in  reasonable  agreement  with  the  more 
accurate  numerical  calculations  of  Sheen  et  al.  [4]  for  stripline,  so  we  will  use  the  analytical  expressions 
to  calculate  IMPs.  It  will  be  assumed  that  this  field  distribution  is  not  changed  significantly  by  the 
presence  of  the  junction  or  the  IMPs,  and  that  the  current  is  uniformly  distributed  through  the  film 
thickness.  The  former  is  approximately  trae  for  small  relative  levels  of  intermodulation  power,  and  the 
latter  will  be  approximately  true  for  films  less  than  2A  thick,  such  as  the  films  considered  here. 

The  expressions  of  [6]  yield  the  factor  in  Equation  (7.6)  that  we  seek: 


iw^jLLr^ifr' , 

Vj  Ttbw  V  M'  7  2  2b 


(7.8a) 


Ji(x)^  1.165Ij 
Vi  ~  TtXiwbff^ 


exp 


w  A?'  .  w 

- <x<—  , 

2  2b  2 


(7.8b) 


where  X  is  the  superconductor  penetration  depth. 

Equations  (7.2)  are  easily  generalized  to  the  situation  in  which  the  junctions  are  distributed  along 
z  =  0,  as  they  are  in  the  resonator  in  Figure  7-1.  Since  I  and  flow  through  the  same  area,  2//7^  can 
be  replaced  with  2J(I,x)/J(^  in  Equation  (7.2). 

The  modified  Equations  (7.2)  along  with  (7.8)  supply  all  the  information  needed  to  evaluate  Equa¬ 
tion  (7.6)  for  the  IMS  current  amplitude  /j  arising  from  the  two  fundamental  signals,  each  with  current 
amplitude  I.  The  integrals  over  x,  although  somewhat  unwieldy,  can  be  solved  in  closed  form. 

Several  resonator  stmctures  like  that  shown  in  Figure  7-1  were  fabricated  in  YBa2Cu307_^  (YBCO) 
from  Josephson  superconductor-normal-superconductor  edge  junctions  made  by  Koren  et  al.  [7].  Mea¬ 
surements  of  the  loss  arising  from  the  junction  in  one  particular  stmcture  were  recently  reported  by  Oates 
et  al.  [8].  They  showed  that  the  junction  dominates  the  resonator  nonlinear  behavior  until  current  densities 
exceed  its  by  an  order  of  magnitude.  We  have  made  intermodulation  measurements  on  the  same 
junction  as  that  measured  by  Oates  et  al.,  so  a  direct  comparison  can  be  made  between  our  derived  values 
and  those  in  [8]. 

Figure  7-3  shows  the  calculated  and  experimental  IM3  currents  for  this  resonator  at  a  temperature 
of  71  K.  Power  measurements  were  converted  to  currents  in  Figure  7-3  using  expressions  given  in  [9]. 
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FUNDAMENTAL  CURRENT  (A) 


Figure  7-3.  Fundamental  and  third-order  intermodulation  currents  in  stripline  resonator  at  71  K.  The  theoretical 
curve  has  no  adjustable  parameters. 


We  have  used  the  values  of  R  and  7^  derived  in  [8]  together  with  the  calculations  outlined  above,  so  there 
are  no  adjustable  parameters  in  the  calculation  of  these  IM3  currents.  The  agreement  is  excellent  at  this 
temperature,  suggesting  that  the  approximations  made  in  the  calculation  are  reasonable  and  that  the  RSJ 
model  can  give  accurate  results  for  YBCO  Josephson  junctions. 

In  future  work  we  expect  to  compare  the  loss  and  intermodulation  results  for  these  Josephson 
junctions  in  resonators  at  different  temperatures  and  external  magnetic  fields  to  further  illuminate  the 
physics  of  these  YBCO  junctions.  In  addition,  we  will  determine  whether  the  coupled-grain  model  can 
explain  both  the  loss  and  intermodulation  measurements  from  a  YBCO  film  without  intentional  junctions. 

T.  C.  L.  G.  Sollner 
J.  P.  Sage 
D.  E.  Oates 
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8.  DIGITAL  INTEGRATED  CIRCUITS 


8.1  ULTRAfflGH  DENSITY  PROM 

A  highly  scaled  programmable  read-only  memory  (PROM)  is  being  developed  by  taking  advantage 
of  ultimate  lithographic  resolution.  This  dense  PROM  device  will  be  a  nonvolatile,  once-programmable, 
fast  alternative  to  mass  storage  media  such  as  compact  disk  read-only  memory  (CDROM),  with  zero 
standby  power  dissipation  and  no  moving  parts.  The  memory  structure  incorporates  a  vertically  self- 
aligned  diode  and  voltage  programmable  link  (VPL)  for  a  minimal  memory  cell  area  architecture  just  one 
lithographic  line  and  space  in  each  dimension.  The  densest  arrays  yet  built,  using  conventional  i-line 
lithography,  have  0.6-/xm  lines  and  spaces,  corresponding  to  a  memory  density  of  70  Mbits/cm^,  while 
individual  bits  of  0.25-jUm  diameter  have  been  built  using  phase-shift  i-line  exposure.  The  simple  nature 
of  the  array  makes  phase  shift  lithography  feasible,  so  with  Lincoln  Laboratory’s  unique  193-nm  stepper, 
a  memory  array  pitch  of  0.2  /tm  is  feasible,  corresponding  to  the  CDROM  capacity  of  640  Mbytes  in 
2  cm^  of  silicon.  A  transistor  version  of  this  device  has  been  previously  reported  in  [1],[2].  The  diode 
form  is  a  simpler,  potentially  denser  version. 

The  basic  memory  element  is  a  VPL  and  a  diode  in  series.  Operation  of  the  single  bit  can  be  seen 
in  Figure  8-1.  This  figure  shows  two  traces,  the  first  being  before  programming.  For  all  biases  less  than 
+6  V,  negligible  current  flows.  Above  6-V  forward  bias,  current  starts  to  flow  through  the  link  insulator 
until  at  10-V  bias  in  Figure  8-1  the  insulating  dielectric  breaks  down,  and  a  conducting  filament  forms 
resulting  in  a  permanent  electrical  connection.  The  second  trace  shows  the  current  after  the  link  is  formed, 
which  is  the  characteristic  of  the  diode,  including  reverse  breakdown  at  -13  V  in  this  case.  A  read  cycle 
would  occur  at  a  word  line  bias  of  <  5  V,  so  the  current  difference  between  a  programmed  and 
unprogranuned  bit  is  many  orders  of  magnitude.  The  purpose  of  the  diode,  unnecessary  for  a  single  bit, 
is  to  provide  isolation  in  a  memory  array. 
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WORD  LINE  VOLTAGE  (V) 

Figure  8-1.  Conduction  characteristics  of  a  single  bit,  during  and  after  writing.  The  first  sweep  shows  picoamps 
for  reverse  bias  and  up  to  6-V  forward  bias,  with  writing  at  +70  V.  The  second  sweep  shows  diode  reverse 
breakdown  at  ~13  V  and  forward  current  limited  at  1  mA. 
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In  an  array,  illustrated  in  Figure  8-2,  the  read  and  write  cycles  are  the  same  as  a  single  bit  with 
additional  provisions  for  unselected  lines.  For  the  bit  being  accessed  the  word  line  is  set  high,  forward 
biasing  the  diode,  and  the  bit  line  is  grounded.  For  all  other  bits,  either  zero  or  reverse  bias  is  needed. 
All  other  bit  lines  are  set  high  and  all  other  word  lines  are  grounded.  In  an  n  x  n  array,  this  results  in 
one  bit  forward  biased,  2(n-l)  half-selected  bits  at  zero  bias,  and  the  remaining  nearly  rP-  unselected  bits 
reverse  biased.  The  reverse  leakage  of  a  single  diode  is  negligible,  but  the  sum  of  all  programmed, 
unselected  bits  could  represent  significant  power  dissipation  during  read.  This  current  is  minimized  by 
small  diode  area,  low  sidewall  leakage  due  to  proper  fabrication  techniques,  and  read  voltage  much  less 
than  write  voltage. 

The  starting  material  for  fabricating  this  PROM  is  a  two-layer  epitaxial  n-type  growth  on  p-type 
silicon.  The  diode  is  formed  in  the  top,  lightly  doped  «-type  layer,  and  the  bottom,  heavily  doped  n+ 
region  forms  the  bit  line.  The  first  step  in  PROM  fabrication  is  to  define  the  diode  pillars  (see  Fig¬ 
ure  8-2).  The  current  process  uses  i-line  lithography  and  reactive  ion  etching  in  He-Cl2.  The  minimum 
feature  size  in  arrays  built  thus  far  is  0.6  |tm,  and  the  etch  depth  is  1.0  /tm.  A  second  lithography  step 
and  1.0-jum  silicon  etch  define  the  rP  bit  lines.  The  severe  topography  generated  by  the  two 
silicon  etches  is  then  planarized  using  low-pressure  chemical  vapor  deposition  silicon  dioxide  from  a 
tetraethylorthosilicate  (TEOS)  source  followed  by  chemical-mechanical  polishing  (CMP)  as  described  in 
[1].  The  oxide  layer  is  conformally  deposited  to  fill  the  area  between  the  pillars,  and  the  wafer  is  then 
polished  to  remove  the  oxide  from  the  pillar  tops.  This  planarizing  step  results  in  a  smooth  surface  with 
only  the  tops  of  silicon  pillars  exposed.  An  unmasked  shallow  p*  implant  forms  the  diode  anode  at  the 
top  of  these  exposed  pillars.  The  VPL  insulator  trilayer,  5  nm  each  of  silicon  nitride,  oxide,  and  nitride, 
is  deposited  by  plasma-enhanced  chemical  vapor  deposition,  and  contact  holes  are  cut  through  the  insu¬ 
lators  to  connect  to  the  n"*"  bit  lines  at  the  array  periphery.  Metal  word  lines  and  other  interconnect  are 
then  sputter  deposited  and  patterned. 
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Figure  8-2.  Schematic  view  of  memory  array.  The  space  between  silicon  pillars  is  filled  with  silicon  dioxide,  which 
is  planarized  to  provide  a  fiat  surface  on  which  the  link  insulator  and  metal  word  lines  lie. 
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Some  of  the  critical  steps  include  alignment  of  the  two  silicon  etches,  silicon  pillar  sidewall 
passivation,  and  CMP  endpoint  control.  Misalignment  during  the  bit  line  patterning  leads  to  etching  part 
of  the  previously  defined  pillar.  A  self-aligned  process  is  being  explored  utilizing  two  crossed  grating 
etches  to  define  the  pillars  and  bit  lines.  A  second  critical  processing  step  is  etch  damage  removal  and 
passivation  of  the  silicon  sidewalls,  especially  in  the  region  of  the  pn  junction  between  the  bit  line  and 
substrate.  Without  proper  passivation  there  is  current  leakage  from  the  ri^  bit  lines  to  the  p-type  substrate, 
leading  to  signal  loss  and  power  consumption.  This  leakage  has  been  minimized  by  growing  a  sacrificial 
oxide  to  consume  the  damaged  silicon  and  then  regrowing  a  high-quality  thermal  oxide  to  complete 
passivation.  The  third  critical  processing  area  is  planarization  of  the  TEOS  oxide.  The  silicon  pillar  tops 
must  be  exposed  for  electrical  conductivity  but  must  not  be  eroded  to  maintain  sufficient  diode  reverse 
breakdown  voltage  for  writing.  Both  the  use  of  a  hard  silicon  nitride  layer  to  stop  the  CMP  and  a  post- 
CMP  highly  selective  wet-chemical  etch  are  being  explored  for  endpoint  control. 

Single  bits  at  0.25-;Um  diameter  and  memory  arrays  at  0.6  pxa.  have  been  successfully  fabricated 
and  tested.  Checkerboards  and  other  patterns  have  been  programmed  in  small  arrays.  Programming  can 
be  accomplished  one  word  at  a  time  with  single  12-V,  1-ps  pulses. 

A  new  reticle  set  has  been  designed  that  includes  arrays  with  selection  circuitiy  and  deep-ultraviolet 
lithography  capability.  The  selection  circuitry  will  allow  on-chip  decoding  of  a  4-kbit  memory  array.  The 
bit  selection  is  performed  by  p-channel  MOS  pass  transistors.  These  transistors  are  fabricated  in  the  n 
epitaxial  layer,  and  transistor  isolation  is  provided  by  the  silicon  etches  required  for  the  memory  array. 
The  reticle  set  includes  arrays  with  feature  sizes  down  to  0.3  ^m  using  248-nm  lithography. 

S.  P.  Doran  P.  W.  Wyatt 

H.  I.  Liu  T.  M.  Lyszczarz 
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